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PROGRESS  REPORT 
OF  CONTRACT  NONR  3173 
GEM  CONTROL  SYSTEMS  STUDY 

1,0  INTRODUCTION 

This  report  will  include  a  summarization  of  that  work 
performed  under  Navy  Contract  NONR  3173  from  November  1961,  to 
November  1962, 

The  AiResearch  GEM  Control  Systems  Study  as  it  existed 

early  in-l962-  (as  report-ed  -in  AP-§047-R,  Reference  1-) . 

consisted  of  an  analytical  formulation  upon  which  one  could 
accurately  predict  the  natural  characteristics  of  virtually 
any  vehicle.  This  formulation  was  incorporated  into  a 
digital  computer  program  used  to  calculate  the  steady-state 
characteristics  and  an  analog  computer' simulation  that 
accurately  represented  the  vehicle  in  the  transient  modes 
of  pitch,  roll,  and  heave. 

It  could  be  considered • that  the  primary  aim  of  the  work 
performed  since  the  last  progress  report  has  been  to  establish 
a  complete  analytical  formulation  of  a  typical  ground  effect 
machine  that  would  be  fully  controllable  in  6  degrees  of  freedom. 

The  initial  approach  was  of  course  completely  general  in 
that  it  could  apply  to  virtually  any  peripheral  jet,  compartmented- 
base  vehicle.  Upon  further  study,  however,  it  was  found  that 
the  analysis  was  growing  far  too  complex  for  any  practical 
application.  However,  defining  a  reference  design  vehicle  and 
making  certain  valid  assumptions  allowed  the  formulation  to  be 
radically  simplified. 

,  AP-5061-R  c 
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In  summary,  the  efforts  have  been  extended  toward  the 
accomplishment  of  the  following  tasks s 

(1)  Simplification  and  linearization  of  the  basic  equations 
of  pitch,  roll,  and  heave  motion.  Previously,  there 
existed  in  this  formulation  nonlinearities  that  were 
shown  to  have  very  little  effect  upon  the  response  of 
the  vehicle.  These  nonlinearities  arose  primarily 
from  the  equations  describing  the  transient  behavior 

of  the  base  pressure.  It  is  now  assumed  that  the  rate 
of  change  of  base  pressure  is  zero — i.e.,  the  new 
pressure  is  reached  instantaneously  following  a  per¬ 
turbation. 

(2)  Development  of  the  remaining  three  equations  of  motion 
for  yaw,  "foirward  flight,"  and  "side  slip"  to  complete 
the  6-degree-of- freedom  simulation. 

(3)  Analytically  describe  the  forces  and  moments  that 
partially  comprise  the  above-mentioned  equations  of 
motion  and  the  former  3  degrees  of  freedom,  arising 
from  aerodynamics  due  to  relative  wind. 

(4)  Development  of  an  analytical  representation  of  the 
water  dynamics  beneath  a  vehicle,  including  the  all- 
important  "hump"  speed. 

(5)  Physically  define  a  complete  vehicle,  including  the 

f 

power  source,  the  air-moving  equipment,  and  the 
propulsive  means. 
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(6)  Analytically  describe  the  dynamics  of  the  engine, 
fan-duct  system,  and  propulsion  propellers. 

(7)  A  preliminary  design  and  simulation  of  a  suitable 
control  system.  Including  the  characteristics  of  the 
system. 

(8)  Construction  of  pilots'  controls  that  will  provide 
an  electrical  signal  to  the  analog  computer 
following  a  control  command  input.  This  system 
includes  a  pilot  display  of  pitch  angle,  roll  angle, 

yaw  angle  (relative  heading),  yaw  rate,  forward  velocity, 
side  velocity,  and  an  x-y  plot  of  position. 

(9)  Generate  and  perform  a  test  program  in  order  to 
evaluate  the  over-all  characteristics  of  the  refer¬ 
ence  design  vehicle. 
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2.0  ANALYSIS  FORMULATION 

The  purpose  of  this  section  will  be  to  discuss  the  steps 
leading  to  the  analytical  formulation  and  analog  computer 
simulation  of  the  Reference  Design  Vehicle  that  is  controllable 
in  6  degrees  of  freedom.  The  appendices  at  the  end  of  this 
report  contain  complete  derivations  of  the  analytics  involved 
in  the  final  formulation;  therefore^  no  derivations  will  be 
presented  here. 

2.1  Linearization  and  Simplification 

In  order  to  utilize  the  analytics  as  they  existed  at  the 
beginning  of  this  phase  of  the  study^  it  was  first  necessary 
to  make  certain  assumptions  in  order  to  reduce  the  complexity 
of  the  over -all  formulation. 

2.1.1  Basic  Simplifying  Assumptions 

(a)  The  base  area  is  divided  into  four  equal  areas  such 
that  =  Ag/4,  The  four  areas  or  compartments  are 

defined  by  the  stabilizing  (compartmentation)  jets 
and  the  peripheral  jets.  This  makes  possible  the 
further  assumption  that  at  steady-state  hover  condi¬ 
tions,  the  pressures  (P^^)  in  all  four  compartments 
(i)  are  equal. 
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(b)  The  above-mentioned  base  compartment  pressure  is  an 
average  steady-state  value  (P^)  plus  a  variation  (pj^). 

The  average  steady-state  value  will  be  seen  to  be  less 
than  the  quantity  Wc/Ag  by  an  amount  that  is  equivalent 
to  the  vertical  thrust  produced  by  the  jets, 

(c)  The  fan  air  flow  is  a  steady-state  value  plus  a 

variation 

(d)  The  products  of  variable  terms  in  the  jet  flow  equations 
are  neglected  in  the  peripheral  jet  equations  but  are 
retained  in  the  compartment  jet  equations  because  of 
relative  magnitude, 

(e)  The  effects  of  pitching;,  rolling,  and  heaving  motions 
upon  the  compartment  jet  heights  are  neglected  (h^j^  =  1). 

(f)  The  rate  of  change  of  pressure  in  a  compartment  is  zero-- 
that  is,  the  new  pressure  level  is  reached  instantaneously 
following  a  perturbation.  Previously,  the  simulation 
contained  the  equations  for  the  rate  of  change  of  pressure 
(Pj^),  These  equations  had  to  be  integrated  in  order  to 
reach  the  new  pressure  level.  The  rate  of  change  of 
pressure  was  so  great  that  the  actual  integration  could 
not  be  done  rapidly  enough.  For  this  reason  the  former 
simulation  was  not  in  real  time  but  had  been  "slowed 
down"  by  a  factor  of  lo.  With  the  assumption  that 
pressure  reaches  its  new  value  instantaneously,  the 
simulation  is  now  in  real  time. 
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(g)  The  peripheral  jets  are  assumed  to  be  always  underfed. 

(h)  The  compartment  jets  are  assumed  to  be  always  overfed. 

Equations  (A-38),  (A-39)j  and  (A-40)  are  the  general  equations 
of®motion  for  heave,  roll,  and  pitchy  It  is  noted  that  they  have 
not  only  been  simplified  using  the  above  assumptions  but  also  that 
they  have  been  normalized  and  linearized  to  a  certain  degree. 

At  this  point  the  base  configuration  of  the  reference  design 
was  introduced,  and. the  equations  for  the  rate  of  change  of. 

o 

pressure  (p^),  flow  and  and  body  motions  were  written 

for  this  particular  model. 

In  Section  3  of  Appendix  A  the  completely  linear  equations 
of  motion  in  pitch,  roll,  and  heave  are  developed.  It  is  noted 
that  Equations  (A-92),  (A-96),  and  (A-lOO)  present  the  variations 
in  heave,  pitch,  and  roll  due  to  fan  flow  variations  and  surface 
disturbances. 

2.2  Six-Degree-of-Freedom  Equations 

Appendix  B  contains  the  development  of  the  equations  of 
motion  in  6  degrees  of  freedom.  Section  1  of  Appendix  B  shows  .. 
the  general  equations  of  motion  of  a  body  free  to  translate  in 
three  dimensions  and  rotate  about  3  body-fixed  axes.  In  order 
to  simplify  these  equations,  it  is  assumed  that  the  origin  of 
the  body-fixed  axis  system  coincides  with  the  center  of  mass 
and  also  that  the  body  fixed  axes  are  the  principal  axes  (those 
axes  about  which  the  products  of  inertia  go  to  zero). 
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The  primary  concern  of  Section  2  is  the  analytical  description 
of  the  external  forces  and  moments  that  act  on  the  vehicle  in 
"flight,"  These  forces  and  moments  will  arise  from  the  following: 

(1)  Propeller  Thrust  -  The  reference  design  vehicle  has 

as  its  primary  source  of  propulsion  two  variable-pitch 
propellers  mounted  on  the  "stern"  of  the  craft.  These 
propellers,  when  manipulated  together,  provide  thrust 
for  forward  flight  and  braking.  When  operated  differen¬ 
tially  they  provide  a  moment  about  the  vertical  axis  for 
yawing.  The  propellers  also  produce  a  pitching  nooienC 
about  the  lateral  (y)  axis, 

(2)  Pitching  and  Rolling.  -  These  motions  about  the  longi¬ 
tudinal  and  lateral  axes  produce  forward  and  side  forces 
resulting  from  components  of  lift.  The  force  in  the 
lateral  direction  due  to  rolling  the  craft  will  be  the 
only  force  available  in  this  direction  other  than  aero¬ 
dynamic  forces » 

(3)  External  Aerodynamics  -  The  vehicle  is  subjected  to  a 
number  of  forces  and  moments  arising  from  external 
aerodynamics.  These  aerodynamic  forces  become  appreci¬ 
able  only  when  the  vehicle  travels  at  high  speed.  They 
are,  however,  analytically  described  over  the  full 
range  of  speeds.  They  are  shown  to  affect  the  following 
areas : 

(a)  Lift  -  at  forward  speed  the  cushion  is  augmented 
due  to  an  aerodynamic  lift. 
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(b)  This  lift  acts  forward  of  the  craft  center  of 
gravity,  therefore  producing  a  pitching  moment „ 

(c)  Stability  in  the  yaw  mode  is  achieved  by  aero¬ 
dynamic  forces  acting  aft  of  the  center  of  gravity, 
thus  producing  yaw  moments  that  tend  to  restore  the 
craft  to  a  longitudinal  translation. 

(d)  Parasite  drag  or  air  friction  produces  a  force 
that  resists  motion. 

(e)  An  induced  pressure  occurs  around  the  periphery  of 
the  vehicle  which  can  be  considered  as  a  variation  in 
the  "ambient"  pressure  which  will,  in  turn,  produce 
slight  variations  in  base  pressure. 


(f)  Momentum  drag  results  from  a  combination  of  forward 
vel  city  and  the  air  flow  to  the  fans.  At  forward 
or  side  flight  the  air  mass  that  is  directed  into 
the  fan  intakes  is  forced  to  change  direction,  which 
produces  a  rather  large  drag  force  on  the  craft. 


As  the  exact  shape  and  configuration  of  the  vehicle  are 
not  established,  it  is  pointed  out  that  many  of  the 
parameters  such  as  lift  and  drag  coefficients,  thrust 
characteristics  of  the  propellers,  and  Internal  aero- 
d3niamlcs  are  only  estimates  based  upon  preliminary 
calculations  by  the  Aerodynamics  Group. 


AP-506I-R 
Page  8 


FOHM  27t7A.t 


•  ' 


^iPesearch  Manu^uturing  Division 

Pk0*nix,  Arixetta 


(4)  Wind  -  To  study  the  handling  characteristics  of  the 
vehicle  in  a  cross  wind,  the  term  "relative  wind"  is 
seen  to  include  both  the  effects  of  flight  and  wind. 

The  final  six  equations  of  motion  are  full  range  equations 
such  that  with  the  proper  input  commands  the  vehicle  can  be 
maneuvered  throughout  a  full  range  of  operations.  It  is  noted 
that  these  equations  no  longer  contain  the  terms  describing  an 
ocean  wave  moving  beneath  the  craft.  The  effects  of  these  waves 
on  a  vehicle  were  described  in  the  April,  I962,  progress  report 
(AP-5047-R).  In  the  present  simulation,  however,  it  was  impossible 
to  include  these  terms,  due  to  a  shortage  of  electronic  multipliers 
in  the  analog  computer, 

2.3  Surface  Deflections.  Base  Pressure  Variations  and  Effects  of  , 

Water  Dynamics  upon  Body  Motions  in  the  Transient 

2.3,1  The  simulation  of  a  vehicle  capable  of  maneuvering  over 
water  would  be  Incomplete  without  a  description  of  the  water 
surface  dynamics  due  to  the  presence  of  a  pressure  area  (vehicle's 
cushion  pressure)  over  the  surface, 

J,  Jo  Stoker*  has  presented  the  general  solution  of  the 
deflection  of  the  water  surface  in  a  free  stream  flowing  beneath 
a  pressure  area  of  a  finite  length.  Stoker's  solution  is  in  the 
form  of  an  Integral  of  an  analytic  function  as  presented  in 
Appendix  C. 


*  Reference  3 
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In  order  to  determine  the  surface  displacement  at  discrete 
points  over  the  length  of  the  pressure  area,  it  was  necessary  to 
integrate  the  functiono  However,  integration  of  the  function 
yielded  not  elementary  functions  but  rather  an  infinite  power 
series  upon  successive  integration  by  parts.  For  this  reason 
it  was  necessary  to  perform  a  numerical  integration.  After 
defining  the  integrals  and  their  paths,  Simponson's  rule  was 
applied  to  the  functions  which,  with  the  aid  of  a  digital  com¬ 
puter  program,  yielded  a  set  of  surface  shapes  that  were  a  function 
of  the  free-stream  velocity  and  applied  pressure. 


This  method  of  solution  was,  of  course,  far  too  complex  and 
lengthy  to  be  of  any  practical  use  other  than  as  a  comparison. 
Therefore,  an  equation  was  Introduced  that  very  closely  approxi¬ 
mated  those  surface  shapes  calculated  by  the  computer,  using  the 
exact  solution.  This  equation  can  be  thought  of  as  a  traveling 
wave  of  fixed  shape  moving  with  the  craft.  From  this  equation  two 
expressions  were  obtained.  One  gives  the  mean  vertical  displace¬ 
ment  as  a  function  of  forward  velocity,  while  the  other  gives  the 
mean  slope  of  the  water  surface  as  a  function  of  forward  speed. 

The  water  slope,  as  it  may  be  considered  a  disturbance,  will 
be  included  in  the  pitch  equation.  This  is  the  term  that  analyti¬ 
cally  describes  the  "hump"  that  is  developed  beneath  and  in  front 
of  a  vehicle  as  it  accelerates  from  hover  to  cruise  speed. 

The  second  of  the  expressions,  that  gives  the  mean  displace¬ 
ment  of  the  surface,  will  be  used  only  as  a  descriptive  term  in  a 
frequency-modifying  term. 
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2.3.2  Relative  wind;,  the  combination  of  vehicle  and  wind 
velocities  relative  to  the  surface  produces  induced  pressures 
around  the  periphery  of  a  vehicle.  For  example,  while  the  vehicle 
is  traveling  forward  at  high  speeds,  a  high-pressure  area  is  pro¬ 
duced  at  the  front  of  the  vehicle  due  to  the  velocity  "head"  and 

a  low-pressure  area  is  created  around  the  sides  and  rear  of  the 
craft.  This  effect,  or  the  variation  in  ambient  pressure  around 
the  jets  causes  variations  in  base  pressure. 

Section  2,0  of  Appendix  C  presents  the  expressions  that 
analytically  show  the  effect  upon  body  motions  of  these  varying 
ambients, 

2.3.3  As  a  vehicle  pitches,  rolls,  or  heaves  while  hovering  or 
traveling  at  low  speed,  the  surface  of  the  water  tends  to  "follow" 
the  motion  of  the  vehicle.  This  is  caused,  for  example,  when  the 
nose  pitches  down  and  the  pressure  beneath  the  nose  momentarily 
rises,  causing  an  increased  depression  of  the  water  surface.  In 
effect,  this  produces  a  decrease  in  the  "cushion  spring  rate" 
resulting  in  a  lower  natural  frequency. 

Section  3.0  in  Appendix  G  presents  the  derivation  of  this 
frequency-modifying  term  for  the  three  motion  modes  of  pitch, 
roll,  and  heave.  It  should  be  noted  that  as  forward  speed  increases 
the  effect  of  water  dynamics  upon  body  motion  diminishes.  The  pre¬ 
viously  mentioned  term  for  the  mean  displacement  is  seen  to  be  a 
function  of  forward  speed,  and  it  is  further  noted  that  this  term 
approaches  zero  as  speed  increases.  Therefore,  this  term  is  used 
to  make  the  frequency-modifying  term  a  function  of  forward  velocity 
as  seen  In  Section  3.0  of  Appendix  C. 


AP-5061-R 
Page  11 


^iPeseardi  Manufacturing  Division 


Photnix,  Aritona 


2 . 4  Reference  Design  Vehicle 

As  mentioned  previously^  the  complete  simplification, 
linearization,  and  normalization  of  the  formulation  partially 
depended  upon  defining  a  reference  design  vehicle.  This  section 
will  present  the  following  information  regarding  the  reference 
design  vehicle: 

(a)  Physical  description,  including  shape  and 
general  configuration,  (Figure  2-1) 

(b)  Air-moving  equipment,  including  power  supply, 
cushion,  propulsive  fans,  and  propellers, 

(c)  Control  system  philosophy, 

(d)  Summary  of  the  parameters  and  functions  used 
in  simulation. 

2,4,1  The  base  planform  is  rectangular  and  symmetrical  about 
the  X  and  y  axes.  The  periphery  is  defined  by  peripheral  jets 
that  are  directed  down  and  inward  at  a  45°  angle  with  the 
vertical.  The  base  is  divided  into  four  equal  areas  by  the 
compartmentation  or  stabilizing  jets  which  are  directed  verti¬ 
cally  downward.  It  is  assumed  that  the  center  of  pressure  of 
each  compartment  is  located  at  its  geometric  center. 

The  jets  themselves  are  divided  into  four  groups,  each 
of  which  defines  a  compartment.  The  air  flow  to  these  jets  is 
supplied  by  one  fan  located  above  the  compartment.  With  this 
scheme  the  flow  and  subsequently  the  pressure  in  any  one  com¬ 
partment  is  controlled  by  varying  the  flow  to  the  fan  itself. 
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The  external  shape  of  the  vehicle  is  streamlined  as  much 
as  possible  with  the  only  major  obstructions  being  the  fan  in¬ 
lets,  cockpit,  and  ducted  propellers  mounted  on  the  stern  above 
the  deck  of  the  vehicleo 

The  power  supply  consists  of  a  pair  of  gas  turbine 
engines  which  drive  two  fans  and  one  d'ucted  propeller  each. 

The  engines  are  basically  constant  speed  with  a  droop-type 
governor.  The  speed  may  be  manually  adjusted  by  means  of  a 
fine  speed  set  point.  This  adjustment  is  provided  so  that  the 
pilot  may  trim  in  roll  and  heave  by  adjusting  the  engine  speed. 

The  "lifting"  fans  are  radial-flow,  constant-speed  fans 
which  discharge  into  ducts  leading  to  the  proper  jets.  Control 
of  the  air  flow  is  achieved  by  varying  inlet  guide  vanes  directly 
upstream  of  the  fans. 

Propulsion  is  obtained  by  variable-pitch,  ducted  propellers 
mounted  on  the  stern  of  the  vehicle.  Each  propeller  is  driven  at 
a  constant  speed  by  a  separate  engine, 

2,4,3  Control  System 

The  basic  control  system  consists  of  the  following  pilot- 
operated  components  i 

(a)  Pitch  Control  -  is  achieved  by  differentially 
controlling  the  pressure  in  the  two  front 
compartments  versus  the  two  rear  compartments. 
Pressure  changes  are  achieve"’  by  varying  the 
air  flow  to  the  compartments.  Variable  pitch 
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inlet  guide  vanes  upstream  of  the  fans  provide 
control  over  the  air  flow  to  the  fans.  The  pilot 
then  can  cause  the  vehicle  to  pitch  by  changing 
the  position  of  the  inlet  guide  vanes.  Pitch  trim 
is  obtained  by  a  separate  fine  adjustment  of  the 
inlet  guide  vanes, 

(b)  Roll  Control  -  is  achieved  in  exactly  the  same 
manner  as  pitch  except^  of  course,  the  two  port  com¬ 
partments  are  controlled  differentially  wjLth  the  two 
starboard  compartments.  Roll  trim  is  achieved  by 
adjustment  of  the  fine  speed  set  point  on  the  engines, 

(c)  Heave  Control  -  It  is  not  anticipated  that  continuous 
control  over  heave  will  be  desirable.  Therefore,  con¬ 
trol  will  be  achieved  by  adjustment  of  the  fine  speed 
set  point  on  the  engines, 

t 

(d)  Fore  and  Aft  Speed  Control  -  is  achieved  by  collec¬ 
tively  changing  the  pitch  on  the  two  ducted  propel¬ 
lers,  A  small  amount  of  additional  fore  and  aft 
"thrust"  is  obtained  by  pitching  the  vehicle  and  thus 
obtaining  a  fore  and  aft  component  of  lift, 

(e)  "Side  Slip"  Control  -  Translation  in  a  side-ways 
direction  is  achieved  by  rolling  the  vehicle  and  thus 
obtaining  the  component  of  lift.  This  is  a  small 
force,  however,  since  high  speeds  in  the  side  direction 
are  neither  anticipated  nor  desirable, 

(f)  Yaw  Control  -  is  achieved  by  differentially 
manipulating  the  pitch  on  the  ducted  propellers. 
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2.4.4  Reference  Design  Parameters  (See  Fig,  2-1) 


Physical  Dimensions 

Ag  =  1593  ft® 

W  =  78^140  lb 

c 

I  =  3,086  X  10®  lb  ft® 

=  8,786  X  10®  lb  ft® 

I  =  10,0  X  10*  lb  ft® 

z 

a  =27  feet 

b  =59  feet 

Zp  =  9  feet 

ap  =  7  feet 

a  =  0.842 

F  =  0,657 


Assumed  Parameters 

cj,  =  0,07 
Cl  =  0,08 

^  percent 

rj^  =  67  percent 

b  =  f«eC 

w 

X  =  14.75  feet 
& 

Ya  =  6,75  feet 
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Operating  Conditions  (Steady  State) 
h  =  1,8  feet 

Pi  =  Pa  =  Pa  =  P4  =  45  Ib/ft^ 

P^^  =  40,11  Ib/ft^ 

Vj  =  196.65  ft/sec  (all  jets) 

p  =  0,002425  lb  secVft^ 

Qfi  =  Qfa  =  ^f3  ==  ^£4  =  ftVsec 

Qi  =  Qs  =  Qa  =  Q4  =  4222,9  ft®/sec 

Qi2  =  Q34  -  294,97  ftVsec 

Qi4  =  Q23  =:  5^3,36  ft^/sec 

Xi  =  0,9156  (page  8)  Appendix  A 

Fi  =  0,9696  (page  5)  Appendix  A 

Fgi  “  (page  4)  Appendix  A' 

®io  (p^ge  Appendix  A 

^KA  ~  0=6563  (page  14  of  AP-5047-R^  Appendix  A) 

Miscellaneous  Parameters 


^Qfi 

dPi 


=  -  0,093. 


^Qfi 

SPi 


9 


ft°/sec 


(The  variation  of  fan  flow  with  base  pressure) 
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Ti 

0.00426  sec 

T2 

0 . 187  sec 

W  b® 

c 

3.87 

W^a® 

c 

2.307 

a 

a+b 

0.3139 

b 

a+b 

0.686 

■ 

1.477 

^  =  10.76  1/sec® 

^4  ^  =  13.79  l/sec* 


TT 

P. 


^  o 

IT 


=  0.8162 


j  §  sinQj  =  0.0982 


0.1633 


Fi  =  0.819 
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\ 

I 


in 

i 


Pitch  Equation  Components  (e) 


“0 

=  5-25  rad/ sec 

^0 

=  0,311 

a? 

=  0,649 

^Qfi 

2  he' 

aQfi 

=  17.89  1/sec^ 

=  21.12  1/sec^ 

^  w 

W* 

=  12.06  1/sec^ 

PgAgb-* 

=  23.78 

0,29a 

Ti  /  _  .  ^ 

=  0«035 


K-  =  0o258  (1  +  e  -  ^  cosa)  = 

y  (X 

(U^)  -  6  (-^  ”  cosa) 
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Fore-Aft  Equation  Coniponents  (x) 
Cd  Pg  Ag  b 


“5r 


-  0.00328  1/ft 


Qp  pg 

-  =  0.01875  l/sec 

c 


I  =  0.458 


I  X.  0.546  l/3ec‘ 


=  0.0333  1/sec^ 

Side  Slip  Equation  Components  (v) 
Cd  pg  Ag  b 


Qp  P8 
17 - 


0.00328  1/ft 


=  0,01875  1/sec 


I  -.2.185 


^  =  0.159  1/sec® 
a® 

Yaw  Equation  Cocnponeots, -(-g/) 


W  g  a 

0<.88l  l/sec' 
z 


^  "  0,00415 
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2.4.5  Control  System  Parameters 

Let:  =  inlet  guide  vane  angle 

0^  =  pitch  command  input  angle 
4)^  =  roll  command  input  angle 
Assumed  Control  System  Functions 

Flow  to  Inlet  Guide  Vane  Angle 

SQ. 


SQi 


=  48 


ft®/ 


sec 


deg  (input] 


Pitch 


^“i  1 

^  ^  ^  +  ^cp 

„  '  ^“i  1  . 

^  ”p 

dtti 

^  IT  +  -^-sT  \ 

dtti  ^ 

^  IT.T-sJ-  ^ 
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%=  In 

=  0«1  sec  (10  rad/sec  lag) 

48 

“  n  +'  6.i»)  ®c 

“tl  +*^.lf)  *c 
48 

'>=’  ■  ri  +  6.U) 

*  11  /lu) 


Roll 

TV.  _  1 

~  +  ^cp 

1 

^  +  -T-cp  s) 


as 


— ^  1 
^  '*■  *  ^cp 


a4  =  — 


1 

^  +  "cp 


^  .=  0.214  deg  (output) 
Cleg  ^ Input) 


'^Q  -  sec  (10  rad/sec  lag) 
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qi 


qa 


qa 


q* 


10,3 

TT+T 

JTTsI 

(1  +  2 

5.1s) 

10.: 

J  4  , 

11  +  1 

)TTsT 

)  ,  . 

a  + 1 

).ls) 

Heave 

It  Is  assumed  that  a  change  In  engine  ay— d'  4t9t  point 
Is  Instantaneous  following  the  pilot  command  (Z  )  <. 

Yaw 

Denote  ^  as  propeller  pitch  angle  and  as 
the  pilot  command 

=  4  deg  (output)  : 

S^cp  deg  I input) 

The  control  Is  assumed  to  have  a  10  rad/sec  lag. 

(ryj.  =  0,1  sec) 

Forward  Velocity 

Denote  as  pilot  command  for  forward  velocity 


=  0.1  sec 
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2 „ 4 „ 6  Engine  and  Fuel  Control 

The  engines  are  shaft,  turbine  engines  with  a  free-type 
power  turbine.  They  are  in  the  1300  to  1600  hp  class. 

In  order  to  eliminate  the  cotnplex  representation  of  such 
an  engine  and  its  fuel  control,  it  was  assumed  that  the 
simplified  system  as  shown  in  Figure  2-2  was  sufficient  for 
this  purpose. 


The  particular  parameters  for  this  system  are  given  below. 


t:. 


0.056 


Ib/hr 

rpm 


SNgp 

W 


12.3 


rpm 

I'b/hr 


_  ')  o447 


=  -OoOOS 


ft/lb 

rpm 


J  ==  0.905  ft  lb  sec^  (Total  inertia  of  rotating 

machinery  referred  to  output  shaft  of  engine) 


■fc  "  -.05  sec 

Tgp  =  1  sec 
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Functions : 

Thrust  (T)  for  a  ducted  propeller  is  a  function  of 
propeller  pitch  angle  forward  velocity  (Oj^)  and 

engine  speed  (N^)o  For  purposes  of  this  simulation  it 

was  assumed  that  T'  vs  could  be  represented  as  shown  in 
Figure  2-3o  For  a  given  propeller  pitch  angle  the  correspond 
ing  T  is  then  obtained  from  the  following; 
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1  .  I'  (1.9  -  0.568  li^)  [i  -11-  (N,  -  N,„)] 

where  is  the  steady  state  operating  speed, 

=  0,485  X  lO  ®  1/rpm 

The  shaft  horsepower  absorbed  in  producing  the  required 
thrust  is  given  in  Figure  2-4, 
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FIGURE  2-3:,  THRUST  VS  PROPELLER  PITCH  ANCLE 
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FIGURE  2-4.  DRIVING  HORSEPOWER  VS  THRUST 
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3o0  DESCRIPTION  OF  GEM  CONTROL  STATION  DESIGN 

The  purpose  of  designing  and  building  a  GEM-control  station 
was  threefold: 

Ic  To  provide  the  physical  means  for  analog  simulation 
of  GEM  handling, 

2.  To  determine  the  specific  performance  characteristics 
of  the  GEM  model  under  study, 

3,  To  provide  a  test  setup  for  selecting  the  best  suited 
control  station  arrangement. 

The  functions  to  be  controlled  by  the  operator  are: 


(a)  Yaw 

(b)  Roll 

(c)  Pitch 

(d)  Fore  and  aft  power 

(e)  Hover  height 


In  order  to  satisfy  point  3,  a  total  of  six  control  units 
have  been  Incorporated  In  the  design.  Figure  3.1  shows  the  control 
station  that  has  been  designed  and  built  to  simulate  the  control 
arrangement  of  a  cockpit. 
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Manual  controls  to  serve  the  control  functions  are  incor¬ 
porated  in  a  pivoted  column  headed  by  a  steering  wheel »  To 
accomplish  the  desired  control  functions  this  device  can  be 
manipulated  by 

(a)  Moving  the  entire  column  fore  and  back  (Fig,  3.4) 

(b)  Swinging  the  wheel  about  the  pivot  axis  of  the 
column  (Figo  3.4) 

(c)  Rotating  the  wheel  about  its  own  axis  (Fig.  3.2) 


fl 


i] 

0 

0 

D 


An  alternative  means  of  steering  is  given  by  a  foot-operated 
treadle  (Fig.  3.5). 


Two  lever  type  positioners  have  been  provided  to  take  care 
of  the  less  frequently  operated  controls  of  the  vehicle  (Fig.  3.2), 

The  mechanical  input  signals  will  be  converted  into  pro¬ 
portional  voltage  signals  by  20  kfi  potentiometers  which  are  for 
trimming  purposes  connected  in  series  with  an  adjustable  potentio¬ 
meter  of  2  kfi. 

The  potentiometers  are  rotated  by  a  pulley  about  which  a  wire 
attached  to  the  control  handle  is  wound.  An  equalizing  extension 
spring  eliminates  any  backlash. 


[ 

|] 

0 

I 

[ 


The  frame  of  the  control  station  has  been  designed  to  accommodate 
any  kind  of  seat--i.e.,  an  ordinary  chair  for  the  operator  of  the  . 
simulator  (Fig.  3.6). 

t 

To  enable  adjustment  of  the  relative  steering  wheel  position, 
the  respective  bars  have  been  hinged  to  permit  relocation  of  the 

t 

steering  wheel  column  (Fig.  3.3). 
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FIGURE  3-1.  GEM  CONTROL  STATION,  GENERAL  VIEW 
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FIGURE  3-2.  gem  CONTROL  STATION,  FRONT  VIEW 
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FIGURE  3-3.  GEM  CONTROL  STATION,  REAR  VIEW 

AP-506I-R 
Page  35 

MF-1197 


'r.OMM  a7t7A.  1 


^iOsseatch  Manufacturing  Division 


•  I 


FIGURE  3-4.  GEM  CONTROL  STATION,  COLUMN  MOVED  BACK  AND  ROTATED 

AP-.506I-R 
Page  36 


FORM  27t7A.  1 


Try 


^iPeseardi  Manufacturing  Division 


Wmwit, 


FIGURE  3-6.  GEM  CONTROL  STATION,  ANALOG  COIPUTER 
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4.0  TEST  RESULTS 

The  six  degrees  of  freedom  analog  computer  simulation  of 
a  completely  controlled  reference  GEM  has  two  purposes.  The 
first  is  to  demonstrate  the  analytical  evaluation  of  GEM's 
generally.  The  second  is  to  show  what  characteristics  will 
be  attained  by  the  reference  design.  The  vehicle  character¬ 
istics  have  been  classified  under  (1)  natural  motions, 

(2)  response  to  commands,  and  (3)  controllability. 

Natural  motions  are  the  transient  motions  of  the  vehicle 
with  locked  controls  following  a  disturbance  such  as  a  sudden 
change  in  a  control  imput.  Response  to  commands  means  how 
much  and  how  fast  a  given  vehicle  attitude  or  velocity  changes 
following  a  control  input  which  commands  this  change.  Con¬ 
trollability  means  how  well  a  human  operator  can  drive  or  pilot 
the  vehicle  through  a  prescribed  or  desired  maneuver. 

The  manner  in  which  specific  data  was  obtained  will  be 
discussed  below  as  the  data  is  presented. 

4.1  Natural  Motions 

The  natural  motions  that  were  recorded  are  heave,  pitch, 
roll,  and  yaw  motions.  The  condition  simulated  was  that  of  a 
vehicle  over  smooth  water  responding  to  step  changes  in  pilot 
commands.  Particular  attention  was  given  to  the  coupling  of 
forward  speed  to  pitch,  side  motion  to  roll,  and  yaw  to  side 
velocity. 
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4 , 1. 1  Heave 


In  the  reference  vehicle,  hovering  height  is  controlled  by 
engine  speed.  This  is  done,  in  preference  to  inlet  guide  vane 
changing,  to  keep  the  engines  at  an  efficient  operating  point. 
This  means  that  the  control  command  for  a  new  hovering  height 
is  a  change  of  engine  speed  set  point. 

Figure  4-1  shows  some  typical  curves  of  engine  speed  and 
vehicle  height  following  a  step  change  in  engine  speed  set 
point.  It  is  seen  that  6  to  8  seconds  is  required  for  the 
height  to  change,  that  there  is  no  noticeable  overshoot,  and 
that  height  follows  engine  speed  very  closely.  The  transfer 
functions  for  the  engine  and  vehicle  show  the  engine  speed 
response  to  set  point  to  have  a  second-order  lag  at  0,64 
radian  per  second  with  a  damping  ratio  of  0.83.  This  is  the 
slow  transient  seen  in  the  traces  (Figure  4-1).  The  vehicle 
response  to  fan  flow  changes  has  a  second-order  lag  at  3  to  4 
radians  per  second  with  a  damping  ratio  of  0.31«  Since  the 
vehicle  is  faster  than  the  engine,  it  would  be  expected  to 
follow  the  engine  speed  rather  closely- -as  seen  in  the  traces. 
Thus,  how  fast  vehicle  height  can  be  changed  is  limited  by  the 
engine  response. 

The  traces  also  show  the  steady-state  gain  variation  with 
forward  velocity.  The  traces  show  the  following  gains  in  feet 
per  lOO-rpm  speed  change; 


Velocity,  fps 

0 

20 

50 

Gain,  ft  per  lOO  rpm 

0.21 

0.24 

0.14 

This  change  is  primarily  due  to  the  change  in  effective  water 
surface  stiffness  with  forward  speed. 
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HOVERING  HEIGHT  VS.  TIME  FOLLOWING  A  SPEED  GQMfAND 


FIGURE  4-1 


AP-5061-R 
Page  41 


|7l7Att 


^iPeseatch  Manu^uring  Division 


4ol.2  Pitch 

The  pilot  controls  pitch  by  means  of  differential  inlet 
guide  vanes.  Figures  4-2(a)j  (b)^  and  (c)  show  pitch  angle, 
vehicle  height  changes,  and  forward  velocity  vs  time  following 
a  pitch  command.  Two  things  can  be  seen  from  these  traces. 

The  first  is  the  frequency  and  damping  of  pitch  (and  heave) 
oscillations.  The  second  is  the  influence  of  pitch  on  forward, 
speed. 


From  the  traces,  the  natural  frequencies  and  damping  ratios 
appear  to  be  the  following: 


Speed 

.  0 

•  20  fps 

4o  fps 

3,93 

4,7 

4,62 

h 

0,42 

- 

- 

The  sloping  curve  at  20  and  4o  fps  as  the  oscillations  damp  out 
precludes  any  reasonable  determination  of  the  damping  ratio. 
However,  it  is  interestii^  to  note  that  the  natural  pitching 
frequency  over  a  hard  surface  with  no  coupling  effects  is  5«25 
radians  per  second  with  a  0,31  damping  ratio.  At  hover.  Figure 
4-2(a),  the  frequency  is  nearer  the  natural  frequency  of  heave 
motions,  which  is  3,95  radians  per  second.  At  20  and  40  fps, 
pitch  still  has  a  strong  effect  on  height  but  the  frequency  is 
nearer  that  of  natural  pitch  oscillations. 

The  forward-speed  traces  show  speed  to  be  affected  by  pitch. 
However,  the  effect  is  seen  primarily  as  the  result  of  pro¬ 
longed  pitching  in  one  direction.  The  oscillations  seen  in  both 
the  pitch  and  the  height  traces  are  virtually  undetectable  in  the 
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speed  traces.  In  Figure  4-2 (a)  the  speed  change  is  almost  linear 
with  time  when  the  pitch  angle  is  constant.  This  is  the  ex¬ 
pected  result  of  constant  acceleration,  since  there  is  virtually 
no  aerodynamic  drag  below  10  fps.  The  slope  of  these  curves 
indicates  an  acceleration  rate  of  about  0„045  g  when  the  pitch 
angle  is  about  3  degrees,  or  0,052  radian.  In  Figure  4-2(b) 
and  (c),  accelerations  are  not  as  readily  seen  due  to  the 
variation  of  slope  with  time.  However,  where  points  were 
checked,  the  slope  changes  were  found  to  be  very  nearly  equal 
to  the  pitch  angle  changes  in  radians. 
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FITCH  AND  FORWARD  SPEED  VS. 
TIME  F(H.LOWING  A  PITCH  COM4AND 

FIGURE  4-2(a) 
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20  FPS  FCXIWARD  SPE^ 

PITCH  AND  FORWARD  SPEED  VS. 
TIME  FOLLOWING  A  PITCH  COMiAND 

FIGURE  4-2(b) 
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40  FPS  FORWARD  SPEED 

PITCH  AND  FORWARD  SPEED  VS. 
TIMS  FOLLOWING  A  PITCH  COMIAMD 

FIGURE  4-2(c) 
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4„1„3  Roll 

Figure  4-3  shows  roll  angle  and  sidewise  displacement  vs 
time  as  roll  oscillations  die  out  naturally.  To  see  the 
oscillation  of  side  motions  due  to  roll  oscillations  it  was 
necessary  to  run  the  traces  simultaneously  on  two  x-y  plotters. 
The  figure  shown  is  a  composite  of  the  two  plots.  The  run  was 
made  by  putting  in  a  full  positive  roll  command  followed  by 
a  full  negative  roll  command  and  rxxtorlng  the  roll  command 
to  zero.  As  seen  in  the  figure,  this  effectively  "rocks  the 
boat"  with  an  amplitude  of  a  few  degrees.  The  frequency  of 
oscillation  determined  from  these  traces  (at  hover)  is  3.3 
radians  per  second  with  0,22  damping  ratio.  Over  a  hard  sur¬ 
face  the  frequency  and  damping  ratio  would  have  been  3  <>93 
radians  per  second  and  0,18,  respectively. 

Figure  4-3  also  shows  that  side  motion  is  oscillatory 
when  the  roll  angle  oscillates.  However,  the  amplitudes  shown 
on  the  traces  are  less  than  0,05  feet  (0,6  inch)  for  the 
indicated  roll  amplitudes  of  about  4  degrees.  As  the  ampli¬ 
tude  decreases  the  side  displacement  lags  the  roll  angle  by 
very  nearly  180°,  This  means  that  during  roll  motions  the 
center  of  rotation  is  not  at  the  center  of  gravity  but  is 
rather  above  it  by  the  ratio  of  side  displacement  amplitude 
to  roll  amplitude.  Although  these  amplitudes  are  not  easily 
obtained  from  damped  transient  data,  the  center  of  rotation 
appears  to  be  about  0,7  foot  above  the  center  of  gravity.  It 
is  expected  that  this  distance  would  be  independent  of  grots 
weight.  Theoretically,  it  should  be  g/co*®  for  harmonic 
rolling  motions.  This  gives  2.93  feet  at  the  natural  fre¬ 
quency  of  3.31  radians  per  second. 
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FIGURE  4-3.  NATURAL  ROLL  AND  SIDE  MOTIONS  VS  TIME 
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The  curves  shown  in  Figure  4-4  show  the  vehicle  in  a 
"fishtailing"  motion.  The  GEM  was  first  broughtup  to  a 
steady  speed  and  then  a  yaw  command  was  given.  After  a 
suitable  yaw  rate  was  obtained#  the  yaw  command  was  reduced 
to  zero.  The  subsequent  yaw  oscillations  are  due  to  aero¬ 
dynamic  moment  resulting  from  side  velocity.  The  curves  show 
that  yaw  causes  a  transfer,  of  velocity  from  to  Uy  and  vice 
versa.  They  also  show  the  period  of  oscillation  to  be  con¬ 
siderably  less  at  100  fps  than  at  50  fps  forward  speed.  At 
50  fps  one  period  of  side  motion  (beginning  at  19  fps)  re¬ 
quired  about  84  seconds.  At  100  fps  one  period  of  side 
motion  (beginning  at  IJ  fps)  required  only  about  58  seconds. 

It  is  seen  also  that  the  period  gets  longer  as  the  amplitude 
of  side  velocity  decreases.  This  is  the  expected  result#  since 
the  aerodynamic  forces  are  proportional  to  the  square  of  the 
velocity. 
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4,2  Response  to  Commands 
4o2„l  Acceleration  and  Braking 


The  data  In  Table  4-1  was  obtained  by  accelerating  and 
braking  the  GEM  with  full  forward  thrust  and  full  reverse 
thrust  commands,  A  stralght-llne  course  was  used. 


Speed  Range 

Acceleration 

Braking 

feet/second 

Time 
(sec, ) 

Dlst, 

(feet) 

^AV 

(g-unlts) 

Time 
(sec, ) 

Dlst. 

(feet) 

^AV 

(g-unlts) 

u 

0  to  17,7 

4,4 

40 

0,125 

18,0 

150 

0.03 

>  53 

0  to  47,5 

16,28 

460 

0,09 

31,8 

600 

0.045 

49,5  to  94,4 

38,6 

3000 

0,035 

19,6 

1400 

0.070 

0  to  4o 

10,8 

250 

0,115 

38,0 

670 

0  to  74,4 

27,0 

1250 

0,085 

61,2 

2000 

TABLE  4-1  ACCELERATION  AND  BRAKING  DATA 

The  first  row  In  Table  4-1  shows  the  low-speed  acceleration 
capability.  The  Indication  Is  a  maximum  accelerating  thrust  of 
0,125  g,  or  12,5  percent  of  the  vehicle  gross  weight.  This 
amounts  to  about  10,000  pounds.  The  maximum  braking  thrust  Is 
0,03  g,  or  2400  pounds.  The  second  row,  which  gives  accislsratlon 
and  braking  through  the  "hump"  speed.  Indicates  less  average 
thrust  for  accelerating  but  more  for  braking.  The  effect  of 
water  surface  response  to  vehicle  motion  Is  seen  In  the  pitch 
attitude.  The  computer  traces  from  which  the  data  In  Row  2  was 
taken  are  reproduced  In  Figure  4-5, 
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Consider  now  the  pitch  angle  and  forward-speed  traces 
in  Figure  4-5.  When  thrust  is  first  applied  (at  zero  forward 
speed)  the  water  surface  is  essentially  flat  and  the  vehicle 
pitches  nose  down  due  to  the  thrust  vector  being  considerably 
above  the  center  of  gravity.  The  nose-down  attitude  adds  to 
the  foirward  thrust.  As  speed  increases,  the  wake  develops  with 
a  wave  crest  forming  at  the  nose  and  a  trough  toward  the  rear 
of  the  GEM.  The  result  is  a  nose-up  pitch  attitude  as  the 
vehicle  follows  the  surface  slope.  The  GEM  begins  to  pitch 
upward  when  the  speed  is  about  16  feet  per  second  and  reaches 
its  highest  nose-up  attitude  (about  3  degrees)  when  the  velocity 
is  about  25  feet  per  second.  This  is  seen  to  be  the  lowest 
slope  or  minimum  acceleration  point  on  the  particular  curves 
shown.  After  thla  the  GEM  gets -over' the  hump  and  begins  to 
run  ahead  of  the  wake  so  that  it  feels  an  increasingly  flat 
surface.  The  pitch  angle  is  essentially  zero  when  the  velocity 
reaches  44  feet  per  second.  At  this  speed  it  is  also  seen  that 
the  slope  begins  to  decrease  as  speed  increases,  showing  that 
the  aerodynamic  drag  begins  to  be  significant.  The  third 
row  in  Table  4-1  shows  that  due  to  aerodynamic  drag  (including 
thrust  loss  due  to  forward  speed)  the  average  accelerating  thrust 
from  50  to  95  feet  per  second  is  only  .0.035  or  about  one 
fourth  of  that  available  at  low  speeds.  Braking,  however,  is 
helped  considerably  at  high  speeds.  It  is  seen  that  the  vehicle 
can  slow  down  twice  as  fast  as  it  can  accelerate  from  50  to  95 
feet  per  second. 
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The  data  in  the  last  two  rows  of  Table  4-1  were  taken 
to  show  the  extent  to  which  water  surface  response  affects 
acceleration  and  braking.  It  is  seen  that  the  average 
acceleration  from  0  to  40  feet  per  second  is  0.09  g  compared 
with  0.115  g  when  the  surface  response  is  included.  The 
figures  in  these  two  rows  may  be  taken  as  "over  land"  values. 
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4.2,2  Side  Motions 


Side  velocity  is  controlled  only  by  rolling  the  GEM. 

The  roll  angle  attainable  is  determined  by  the  maximum  ex¬ 
cursion  of  the  inlet  guide  vanes  and  is  about  0,03  radian 
or  1.72  degrees.  At  this  roll  angle  the  GEM  was  accelerated 
sideways  at  various  forward  velocities.  The  results  are 
summarized  in  Table  4-2 . 


(fps) 

0 

(rad) 

U 

y 

Range 
( fps ). 

Time 
(sec, ) 

^AV 

G-Units 

0 

0,030 

0 

to 

+10 

12,5 

0,025 

0,83 

+10 

to 

-10 

21,8 

0,0285 

0,95 

0 

to 

+20 

31,8 

0,0195 

0,65 

+20 

to 

0 

15o6 

0,040 

1,33 

4o 

0,022 

0 

to 

+10 

18 

0,0173 

0,79 

+20 

to 

0 

25 

0,025 

1,14 

60 

0,019 

0 

to 

+10 

13o5 

0,023 

1.21 

0 

CVJ 

+ 

to 

0 

28,5 

0,022 

1,16 

80 

0,022 

0 

CVl 

+ 

to 

0 

27.5 

0,023 

1.05 

SIDE  MOTIONS 
TABLE  4-2 
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The  data  in  Table  4-2  shows  several  characteristics  of  side 
motions.  The  zero  forward-velocity  data  shows  that  more 
acceleration  is  available  for  reducing  than  for  increasing 
velocities.  This  is  the  expected  result  due  to  aerodynamic 
drag.  It  shows  also  that  acceleration  decreases  as  velocity 
increases.  The  maximum  attainable  side  velocity  is  about  20 
feet  per  second.  This  appears  to  be  about  the  hump  speed. 

However,  it  actually  is  the  result  of  a  balance  between  side 
force  due  to  roll  and  side  aerodynamic  drag,  because  the  water 
surface  response  due  to  side  motions  was  not  included  in  the 
analog  computer  simulation. 

The  variation  of  side  acceleration  with  speed  appears  to 
be  somewhat  erratic.  The  traces  indicate  that  side  velocity 
has  a  somewhat  oscillatory  character  when  forward  velocity  is 
present.  Figure  4-6  has  been  included  to  show  this.  The 
reason  for  this  somewhat  erratic  side  velocity  is  that  the 
center  of  pressure  for  side  drag  is  aft  of  the  center  of  gravity, 
causing  the  craft  to  yaw,  and  yawing  causes  an  additional  side 
acceleration  due  to  coupling  with  forward  velocity.  The  operator 
attempted  to  hold  yaw  rate  zero  during  this,  test  and  succeeded 
in  causing  it  to  oscillate  about  the  zero  yaw  rate  point.  The 
result  is  that  the  side  accelerations  due  to  roll  only  could 
not  be  separated  from  the  coupling  effects.  Therefore,  the  data 
in  Table  4-2  at  forward  velocities  can,  at  best,  show  trends  only. 
It  appears  that  side  acceleration  is  not  much  affected  by  forward 
speed  except  through  yaw. 

Although  pilot  opinions  will  be  discussed  in  Section  6,  it 
is  worth  mentioning  here  that  most  pilots  would  prefer  to  have 
more  side  force  available  for  maneuvering  purposes. 
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4o2.3  Holding  Course  in  a  Wind 

Since  the  GEM  has  no  solid  connection  with  the  surface, 
it  is  very  susceptible  to  being  blown  off  course  by  surface 
winds.  Two  means  are  available  for  coping  with  undesirable 
wind  effects.  The  first  is  to  roll  into  the  wind  so  that  the 
side  force  due  to  base  pressure  equals  the  side  drag.  The 
second  is  to  yaw  the  GEM  into  the  wind  so  that  a  component 
of  the  thrust  vector  balances  the  side  aerodynamic  force. 

The  main  forces  involved  in  rolling  into  the  wind  are 
illustrated  in  Figure  4-7. 


ROLLING  INTO  A  WIND 
FIGURE  4-7 

In  performing  this  maneuver  it  was  found,  as  expected,  tliat  the 
pilot  had  no  particular  difficulty  holding  zero  side  speed  rela¬ 
tive  to  the  ground  in  side  winds  of  10  and  20  feet  per  second 
(fps).  This  test  was  run  at  0,  10  fps,  and  20  fps  forward 
velocities.  In  some  cases  the  pilot  held  zero  side  velocity 
from  the  start;  in  others  the  GEM  was  allowed  to  drift  with  the 
wind  velocity.  It  was  then  rolled  into  the  wind  to  reduce  the 
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velocity  to  zero,  much  as  in  the  side  acceleration  runs  mentioned 
previously.  The  traces,  in  fact,  closely  resemble  those  shown 
in  Figure  4-6  if  they  are  interpreted  in  terms  of  air  speed  rather 
than  ground  speed.  Another  effect  that  should  be  noted  here  is 
that  a  small  differential  thrust  is  needed  to  prevent  yaw  due  to 
the  aerod3mamic  moment.  This  again  is  similar  to  the  "side 
acceleration"  runs  in  that  the  pilot  tends  to  generate  slow 
oscillations  in  finding  the  proper  yaw  command. 


The  main  forces  involved  in  yawing  into  the  wind  are 
Illustrated  in  Figure  4-8. 


a )  VOLOCITIES 


b)  FORCES 


Uy  =  cos  f 


C )  RESULTANT  VOLOCITIES 

YAWING  INTO  A  WIND 
FIGURE  4-8 
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Figure  4-8(b)  shows  the  forces  as  they  are  simulated  on  the 

computer.  Since  the  drag  force  is  colinear  with  the  thrust, 

T  +  these  two  forces  must  be  equal  for  constant  velocity  U  . 
p  s  ^ 

That  is 


If  U  =  const  then  T  +  T  -  F 

X  p  S  X 

Furthermore,  since  there  is  only  one  force  in  the  Y  direction 
it  must  be  zero  for  constant  velocity.  This  means  that 

U„  =  -V„  cosil/  if  U„  =  constant 
y  w  ^  y  . 

The  result  of  these  two  conditions  is  that  the  resultant 
velocity  will  also  be  constant  but  may  not  be  in  the  desired 
direction.  In  any  case, 

Uj^  =  cos^f 

V„  cos^ 

a  =  ^  -  tan  ^  — w - 

X 


The  conditions  that  have  a  prescribed  value  U  and  be  directed 
perpendicular  to  the  wind  are  the  following: 

(a)  a  =  0  or  Ujj  tan^  =  cosf 

(b)  cos^ 

That  is  tan^ 


2  j  Z  V 
COS  t  +  sin  ^  _  w 

sin^V'  ” 
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(c)  Tp  +  Tg  =  (U  cosf  +  sin^/) 

Attempts  to  fly  a  straight  course  in  a  wind  without  rolling 
indicate  that  it  can  be  done  but  that  it  takes  several  minutes 
to  find  the  proper  control  setting  because  of  the  slow  yaw 
transients.  It  was  considered  a  commendable  feat  when  a  pilot 
succeeded  in  holding  a  fixed  position  in  space. by  yawing  the 
craft  90  degrees  into  a  4o-fps  wind.  Subsequent  experience  in 
flying  prescribed  courses  further  pointed  out  the  difficulty 
of  controlling  yawing  position  with  the  long  transient  times 
involved.  It  is  suggested  at  this  point  that  greater  yaw 
stability  be  built  into  future  designs.  Physically,  this 
means  more  vertical  stabilizer  area.  Analytically,  it  means 
moving  the  center  of  pressure  farther  aft.  It  was  located 
6  feet  behind  the  center  of  gravity  (about  10  percent  of  the 
vehicle  length)  in  the  simulation  discussed  here. 

4.2,4  Turning 

Turning  is  done  for  the  purpose  of  altering  the  vehicle 
heading  or  course.  Vehicle  heading  is  controlled  by  yawing. 
However,  unlike  an  automobile,  a  GEM  can  yaw  while  moving  in 
a  given  direction  without  noticeably  changing  its  course.  The 
strong  tendency  to  sideslip  requires  that  side  force  as  well  as 
yaw  moment  be  used  in  turning.  The  yaw  capability  built  into 
the  simulation  was  limited  to  a  maximum  rate  of  lO  degrees  per 
second.  This  is  not  realistic  for  a  real  vehicle  because  the 
limiting  yaw  rate  will  vary  with  forward  speed  due  to  aerodynamic 
changes  and  changing  thrust  gain,  Tha  changing  gains,  however, 
were  realistically  simulated.  The  required  pilot  input  at  zero 
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propeller  thrust  is  high  due  to  low  propeller  gain.  The  required 
pilot  input  at  high  velocities  is  high  due  to  increased  yaw 
stability.  At  moderate  speeds  the  yaw  gain  in  response  to 
pilot  commands  is  maximum. 

Three  turning  situations  were  used  to  demonstrate  turning 
capability.  The  first  was  yawing  in  place  (zero  linear 
velocities).  The  second  was  coordinated  turning  (zero  sideslip). 
The  third  was  a  sideslipping  turn  with  U  approximately  equal 
to  Uy.  Yawing  in  place  is  governed  by  the  built-in  yaw  rate 
limit  of  10  degrees  per  second.  Therefore  the  minimum  time 
to  turn  l8o  degrees  is  18  seconds.  With  a  yaw  rate  meter  the 
pilot  can  readily  hold  any  desired  steady  yaw  rate. 

A  coordinated  turn  is  illustrated  in  Figure  4-9* 


PATH 


FIGURE  4-9.  COORDINATED  TURN 


In  this 
velocity  and 


turn^  Uy  =  0  and  the  yaw  rate  is  related  to  forward 
radius  of  curvature  by 
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force  must  be 
force  Is 


Since  there  is  no  sideslip,  must  be  supplied  by  roll  only. 
From  Section  4.2,2  the  maximum  side  acceleration  is  0,03  g. 
Therefore,  the  minimum  radius  of  curvature  for  coordinated 
turns  should  be 

The  computer  runs  show  a  slightly  smaller  radius  than  this. 
These  are  shown  in  Table  4-3. 


FORWARD 

SPEED 

TURNING  RADIUS 

COMPUTER 

EQUATION 

12  FPS 

145  FT 

149  ft 

20  FPS 

405  FT 

414  FT 

TABLE  4-3.  COORDINATED  TURN  RADIUS 

WITH  MAXIMUM  ROLL  COMMAND 


The  computer  traces  from  which  the  data  in  Table  4-3  were 
taken  are  reproduced  in  Figure  4-10. 


In  order  to  perform  such  a  turn,  a  centripetal 
applied  to  the  vehicle.  The  magnitude  of  this 
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A  sideslipping  turn  Is  Illustrated  In  Figure  4-11. 


PATH 


U  F 


FIGURE  4-11.  SIDESLIPPING  TURN 


In  this  turn,  there  is  side  force  due  to  roll  and  sideslip.  The 
centripetal  force,  which  is  normal  to  the  resultant  velocity  U, 
is  given  by 


Fjj  =  F^  sin  a  +  Fy  Cos  a 


The  tangential  force  is 


F^  =  F^  Cos  a  -  Fy  sin  a 


If  the  magnitudes  of  the  velocities  are  constant,  then  F,j,  must 


be  zero,  giving  F  =  F  tan  a  and 

X  y 


=  ^y 


The  resulting  radius  of  curvature  will  be 

R  U® 
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Numerous  constant  U  and  f2  sideslipping  turns  were  ejceeuted. 

It  was  found  quite  difficult  to  control  sideslip  angle.  Velpcityj 
sideslip  angle>  end  turning  radius  for  several  of  these  turns 
are  listed  in  Table  4-4. 


TANGENTIAL 

VELOCITY 

SIDESLIP 

ANGLE 

TURNING 

RADIUS 

^n/^c 

18,2  fps 

57° 

130  ft 

0,0792 

20.9 

VJl 

0 

167 

O.O8I3 

22.2 

540 

195 

0,0786 

23.3 

55° 

215 

0,0785 

25.3 

51° 

254 

0,0783 

27o3 

49° 

298 

0.0777 

30,3 

46° 

374 

0.0763 

33,6 

45° 

472 

0.0743 

TABLE  4-4.  SIDESLIPPING  TURNS 

In  Table  4-4  it  is  seen  that  a  20.9-fps  turn  can  be  made 
with  a  radius  of  167  feet.  The  minimum  coordinated  turn  radius 
at  20  fps  was  405  feet.  A  sideslipping  turn  with  this  radius 
can  be  made  at  a  speed  in  excess  of  30  feet  per  second.  Table 
4-4  also  shows  an  effective  centripetal  force  equal  to  0.075 
to  0,080  for  sideslipping  turns,  whereas  the  value  for 
coordinated  turns  is  only  0.03  W^.  It  is  apparent  that  over 
half  of  this  is  due  to  the  centripetal  component  of  propulsion. 
Figure  4-12  shows  the  computer  trace  from  which  the  20.9  &nd 
22v2-fps  points  were  taken. 
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4,3  Controllability 

'f 

4,3ol  Introduction 

The  attempt  has  been  made  to  evaluate  seriously  a  vehicle 
being  controlled  by  a  human  operator.  However  excellent  the 
simulation  may  be  In  repxoduclng  or  predicting  natural  vehicle  [ 

motions,  response  to  winds  and  waves,  or  reaction  to  a  control 
command.  It  falls  short  insofar  as  It  falls  to  give  some  answers  to  ( 
the  obvious  question,  "Can  a  man  really  pilot  this  vehicle  Into 
a  harbor  and  dock  It  wltlnout  endangering  other  small  craft  In  * 

the  area?"  If  the  simulation  can  give  answers  to  such  basic 
questions.  It  will  be  of  real  value  In  evaluating  proposed 
vehicles  and  proposed  control  systems.  It  Is  felt  that  the  i 

results  described  In  this  section  show : 

1 

( 

(a)  That  the  present  simulation  gives  considerable  insight 

Into  the  practical  aspects  of  designing  controllable  i 

gem's  and  - 

(b)  That  additionaL  work  is  needed  in  giving  the  pilot  j 

an  immediate  realistic  feeling  for  vehicle  motions 

as  they  occur. 

Several  problems  are  encountered  In  trying  to  put  a  human 
operator  in  command  of  a  vehicle, 

(a)  The  vehicle  must  be  controllable. 

For  the  reference  GEM  an  analysis  was  performed 
to  determine  the  most  feasible  means  of  providing 
pitch  control,  side  force,  propulsion,  braking,  and 
yaw  moments. 
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The  pilot  must  have,  at  his  control  station,  physical 
hardware  with  which  he  commands  the  various  vehicle 
Inputs . 

A  rudimentary  control  station  was  constructed 
as  described  in  Section  3.  It  has  been  called  a 
"universal  control  station"  because  most  of  the 
pilot  motions  that  are  normally  found  in  vehicles 
have  been  included  without  designating  the  particular 
pilot  motion  that  controls  each  vehicle  input. 

The  pilot  must  be  continuously  informed  of  vehicle 
motions,  attitudes,  locations,  etc. 

Initial  trial  runs,  which  were  made  to  check 
out  the  program,  used  a  standard  X-Y  plotter  for 
course  plotting,  a  Sanborn  8-channel  recorder  for 
speeds,  heading,  etc.,  and  one  or  two  dial  meters. 
During  these  initial  runs  it  was  found  difficult  to 
keep  track  of  heading  and  yaw  rate,  so  that  the  GEM 
tended  to  spin  in  rather  tight  loops.  It  was  found 
that  a  yaw  rate  meter  helped  the  operator  very  much. 
The  final  pilot  display  employed  six  meters,  showing 
pitch  angle,  roll  angle,  forward  velocity,  side 
velocity,  heading  (^)  and  yaw  rate,  and  an  X-Y 
plotter  to  show  the  vehicle  path  relative  to  surface 
coordinates . 


rORM  27«7A>1 


4„3o2  Control  Modes 


Basically  the  pilot  controls  four  inputs  to  the  vehicle. 
These  are  thrust  (acceleration  and  braking),  roll  moment 
(producing  roll  angle  and,  thus,  side  force),  pitch  moment 
(producing  pitch  angle  and  thrust),  and  yaw  moment.  With  the 
flexibility  of  a  computer  simulation,  it  was  possible  to  use 
any  of  the  pilot  motions  to  control  a  given  vehicle  motion. 

In  the  actual  runs,  four  control  modes  were  used  to  show  their 
relative  merits  and  the  ability  of  a  pilot  to  learn  new  control 
modes. 

4,3.2,!  Control  Mode  1 

Parameter  Pilot  Motion 


Thrust 

Roll 

Pitch 

Yaw 


Left-hand  hand  lever 
Wheel-base  rotation 
Push-pull  wheel 
Wheel  rotation 


4,3o2,2  Control  Mode  2 


Parameter 


Pilot  Motion 


Thrust 

Roll 

Pitch 

Yaw 


Foot  pedals 
Wheel-base  rotation 
Push-pull  wheel 
Wheel  rotation 
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4.3. 2.3  Control  Mode  3 


Parameter 

Thrust 

Roll 

Pitch 

Yaw 

4. 3. 2. 4  Control  Mode  4 

Parameter 

Thrust 

Roll 

Pitch 

Yaw 


Pilot  Motion 

Push-pull  wheel 
Wheel  rotation 
Righ-hand  hand  lever 
Foot  pedals 


Pilot  Motion 

Left-hand  hand  lever 
Wheel  rotation 
Right-hand  hand  lever 
Foot  pedals 


4. 3. 2. 5  Discussion 

Rotation  of  the  wheel  base  and  push-pull  motion  of  the 
wheel  were  centered  by  springs.  For  this  reason  these  motions 
required  a  continuous  force  by  the  pilot.  The  other  motions 
(wheel  rotation,  hand  levers,  and  pedals)  can  be  moved  to  a 
desired  setting  and  left  unattended  until  a  change  is  desired. 

In  Control  Modes  3  and  4,  the  wheel  base  rotation  was  locked  out. 
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4o3o3  Evaluation  Courses 

Evaluations  of  vehicle  controllability  were  divided  into 
low-speed  and  high-speed  maneuvers.  For  this  purpose  low  speed 
is  defined  as  speeds  below  20  feet  per  second--below  hump  speed. 

High  speed  is  above  hump  speed.  The  high-speed  runs  were 
generally  made  above  40  feet  per  second,  and  speeds  in  the  80- 
to  lOO- feet -per -second  range  were  not  uncommon. 

Three  basic  maneuvers  were  used  to  demonstrate  controllability. 
These  are: 

(a)  Flying  a  straight  course  and  stopping  in  a  prescribed 
area. 

(b)  Flying  a  simple  90-degree  turn  and  stopping  in  a 
prescribed  area. 

(c)  Flying  from  point  A  to  point  B  in  the  presence  of 
obstacles,  including  going  through  a  relatively 
narrow  passage. 

Reduced  reproductions  of  the  four  plates  used  on  the  X-Y  plotter 
are  shown  in  Figures  4-13,  4-14,  4-15,  and  4-16.  The  actual 
size  of  the  plates  used  was  11  by  17  inches,  with  a  10  x  15-inch 
area  of  grid  lines.  The  low-speed  courses  (la,  lb,  and  Ic) 
shown  in  Figures  4-13  and  4-l4  used  a  scale  of  200  feet  to  the 
inch.  The  small  craft  shown  on  these  plates  are  30  Jjy  60  feet, 
which  is  approximately  the  size  of  the  reference  design.  The 
high-speed  courses  shown  in  Figures  4-15  and  4-16  used  a  scale 
of  1000  feet  to  the  inch.  Note  that  the  dotted  rectangular  area 
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at  the  bottom  of  Figure  4-16  is  a  reduced  drawing  of  the  low- 
speed  harbor  course  in  Figure  4-14.  These  two  courses  taken 
together  represent  a  high-speed  approach  to  the  coast  line 
through  a  straight  followed  by  a  low  speed  course  through  a 
breakwater  and  harbor. 
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FIGURE  4-14. LOW  SPEED  HARBOR  MANEUVERING  COURSE  Ic 
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4o3<.3  Results  of  Controllability  Runs 
4 „ 3  o  3  o 1  Procedure 

The  general  procedure  for  flying  the  courses  was  to  allow 
the  pilots  some  practice  time  to  become  familiar,  with  the  control 
mode  and  then  to  fly  all  six  courses  in  that  modeo  Three  pilots 
were  used  throughout  the  test  runs„  Usually,  all  three  pilots 
performed  the  low-speed  runs  before  any  started  the  high-speed 
runSo  In  all  cases,  the  pilot  was  allowed  reruns  at  his  own 
discretiono  It  should  also  be  noted  that  all  of  the  preliminary 
familiarization  with  vehicle  characteristics  was  done  in  Mode  1. 
Each  pilot  had  several  hours  of  flying  time  in  Mode  1  before 
any  runs  were  made  for  the  record „  Following  a  change  in 
control  mode,  the  pilots  generally  made  runs  for  the  record 
beginning  with  their  third  or  fourth  trial  run. 

During  the  test  runs  continuous  recordings  of  forward 
velocity,  side  velocity,  yaw  angle,  roll  angle,  pitch  angle, 
and  hovering  height  variation  were  made.  To  correlate  these 
traces  with  the  X-Y  plotter  course  record,  both  the  Sanborn 
trace  and  the  X-Y  plotter  pen  were  blipped  at  lO-second  intervals. 

The  low-speed  runs  were  made  in  Modes  1  and  3  only.  The 
reason  is  that  all  of  the  pilots  preferred  to  use  pitch  only 
for  low-speed  propulsion.  Therefore,  since  Modes  2  and  4 
differed  only  in  the  thrust  control,  they  were  the  same  as 
Modes  1  and  3  at  low  speeds. 
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4.3o3.2  Flying  a  Straight  Course 


The  tendency  of  the  GEM  to  sideslip,  together  with  the 
slow  yaw  transients,  makes  flying  a  straight  course  more 
difficult  than  It  would  at  first  appear.  Figure  4-17  shows 
the  low-speed  runs  In  Control  Modes  1  and  3  by  each  of  the 
three  pilots.  The  times  required  to  traverse  the  1000  feet 
from  A  to  B  are  shown  In  Table  4-5, 


Pilot  1 

Pilot  2 

Pilot  3 

Mode  1 

66  Sec, 

67  Sec , 

62  Sec, 

Mode  3 

80  Sec,* 

77  Sec, 

150  Sec. 

TABLE  4-5o  Time  To  Fly  Course  la 

As  seen  In  Figure  4-17,  all  three  pilots  were  able  to  hold 
a  straight  course  within  15  or  20  feet,  achieve  an  average 
speed  near  I5  feet  per  second,  and  stop  within  the  vehicle 
dimensions  of  the  preset  point.  This  was  In  Mode  1,  with 
which  all  pilots  had  had  considerable  flying  time.  It  Is 
Interesting  to  note  that  an  Idealized  run  with  acceleration 
at  0,045  g  (the  low-speed  value  due  to  pitch  from  Section  4,1,2), 
constant  speed  20  feet  per  second,  and  deceleration  at  0.045  g 
would  require  64  seconds. 

The  Mode  3  tests,  which  were  run  after  the  Mode  1  tests, 
show  considerably  less  skill  on  the  part  of  the  pilots.  The 
difference  lies  In  the  change  to  using  foot  pedals  for  yaw 
control  and  turning  the  wheel  for  roll  control.  The  first 
run  of  Pilot  1  shows  Initial  inability  to  correct  a  slow 

*Run  2 
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FIGURE  4-17  COURSE  la  TEST  RUNS  -  LOW  SPEED 
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divergence  from  the  prescribed  course,  followed  by  a  futile 
attempt  to  stop  and  recover.  This  is  a  situation  that  is 
neither  uncommon  nor  difficult  in  car  or  truck  driving.  If 
the  driver  finds  he's  missing  the  garage  door,  he  simply 
stops,  backs  up,  and  proceeds  slowly  onto  the  proper  course. 
This  would  be  the  obvious  way  out  for  Pilot  1  in  the  run 
being  considered.  However,  the  display  did  not  give  the  pilot 
an  immediate  feeling  for  course,  heading,  and  the  required 
changes.  This  is  a  pattern  that  recurred  throughout  the  test 
runs.  As  a  rule,  when  a  Pilot  found  his  course  to  be  far 
enough  off  to  require  a  drastic,  immediate  revision,  he 
became  confused  as  to  the  heading  of  his  craft  and  was  unable 
to  achieve  the  required  corrections. 

The  second  run  by  Pilot  1  shows  considerable  improvement 
in  holding  course.  However,  the  time  increase  from  66  to  8o 
seconds  indicate  the  extra  caution  with  an  unfamiliar  control 
mode.  Pilot  2's  run  time  increase  shows  additional  caution, 
while  that  for  Pilot  3  shows  that  he  found  it  very  difficult 
to  make  the  required  corrections  even  though  his  speed  was 
less  than  half  of  that  used  in  the  Mode  1  run  all  three  pilots 
had  difficulty  stopping  in  the  prescribed  area.  The  reason 
for  this  is  a  time  delay.  When  the  Mode  3  runs  were  made, 
several  days  had  elapsed  since  the  previous  low-speed  runs 
and  the  pilots  had  forgotten  the  required  stopping  distance 
(about  140  feet)  which  was  learned  by  experience  during  the 
test  runs. 
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Figure  4-1.8  shows  the  high-speed  runso  The  high-speed 
course  has  a  length  of  7000  feet.,  From  the  acceleration  data 
in  Table  4-1,  an  estimate  can  be  made  of  the  time  to  make  this 
7000-foot  run.  With  an  acceleration  rate  of  0,09  g  from  0  to 
50  fps  and  0,035  g  from  50  to  100  fps,  the  GEM  would  reach 
3760  feet.  Deceleration  at  0,07  g  from  lOO  to  50  fps  and 
0„045  g  from  50  to  0  fps  would  require  57  seconds  and  a 
distance  of  2530  feet.  The  middle  710  feet  would  be  traversed 
in  7  seconds.  This  gives  a  total  estimated  run  time  of  126 
seconds.  Figure  4-19  shows  three  typical  speed  versus  time 
traces  for  the  test  runs  in  Figure  4-18,  It  is  seen  that 
all  pilots  made  their  runs  in  less  than  126  seconds.  The 
initial  accelerations  in  the  three  traces  shown  are  0.136  g, 

0,117  g)  ^tid  0,132  g„  It  is  apparent  that  the  pilots  were 
using  both  full  propeller  thrust  and  nose-down  pitching  for 
maximum  acceleration.  The  braking  curves  likewise  show  that 
pitch  was  used  to  augment  the  propellers.  Counting  the  blips 
in  Figure  4-18  shows  that  most  of  the  runs  were  made  in  about 
loo  seconds. 

The  test  runs  in  Figure  4-18  show  the  results  of  experience. 
In  the  Mode  1  runs,  all  three  pilots  held  good  courses  and 
stopped  very  near  the  center  of  the  target  area  (  a  500-foot 
square).  In  the  Mode  3  runs,  two  pilots  overshot  the  target 
area  by  failure  to  brake  soon  enough.  However,  in  the  Mode 
4  runs,  which  were  very  similar  to  those  of  Mode  3,  all  three 
were  again  able  to  stop  within  the  target  area. 
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4,3.3o3  Performing  a  Simple  Turn 

The  low-speed,  90-degree  turn  tests  by  the  three  pilots 
in  Modes  1  and  3  are  shown  in  Figures  4-20,  -22,  and -24,  The 
Sanborn  traces  corresponding  to  these  tests  are  shown  in 
Figures  4-21,  -23,  and  -25o  On  the  best  run  for  each  pilot, 
small  arrows  have  been  dubbed  in  to  show  the  actual  heading 
of  the  GEM  at  points  along  the  curve. 

In  Figure  4-20,  Mode  1,  the  turn  was  being  made  much  too 
slowly.  Pilot  1  slowed  the  vehicle-“as  shown  between  70  and 
80  seconds  in  Figure  4-21--until  a  90-degree  yaw  angle  was 
reached  and  then  accelerated  toward  the  target  area.  The 
Mode  1  test  shows  an  initial  zero  heading,  with  the  course 
drifting  to  the  left  due  to  sideslip.  The  turn  toward  the 
target  area  is  a  nearly  constant  velocity  sideslipping  turn. 

The  final  10  seconds  of  the  run  shows  an  excessive  yaw  angle 
accompanied  by  sideslip  which  continued  to  carry  the  GEM  away 
from  the  target  area. 

In  Figure  4-22,  Pilot  2  has  a  nearly  perfect  heading  and 

course  initially.  At  40  seconds  he  was  sideslipping  into  the 

turn.  The  roll  angle  between  30  and  4o  seconds  (Figure  4-23, 

Mode  1)  is  the  primary  cause  of  a  turning  course  at  this  point. 

By  50  seconds  Pilot  2  had  yawed  to  a  heading  along  the  courge, 

and  by  60  he  had  developed  considerable  sideslip  to  the  outside 

of  the  turn--which  is  normal  for  a  minimum-radius  sideslipping 

turn.  The  failure  to  continue  on  into  the  target  area  is  a 

combination  of  braking  which  eliminated  the  velocity  component 

in  the  direction  of  the  heading  and  a  reduction  of  roll  angle 

to  zero  which  failed  to  eliminate  side  velocity.  It  is  seen 

in  Figure  4-23  that  U  is  negative  at  70  seconds, 

2% 
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In  Figure  4-24,  Pilot.  3  exhibited  good  control  in  correcting 
an  initial  drift  to  the  righto  After  20  seconds  a  0-degree 
nonsideslipping  course  had  been  achievedo  From  30  to  50 
seconds,  a  small-radius  sideslipping  turn  is  achieved.  From 
50  seconds  onward  he  controlled  the  course  strictly  by  sideslip. 
The  traces  in  Figure  4-25  show  that  the  braking  was  accomplished 
by  a  reduction  of  roll  angle  so  that  sideslip  carried  the  GEM 
into  the  target  area. 

Two  other  observations  can  be  made  regarding  these  tests. 

The  first  is  in  regard  to  turning  radius.  The  smallest  radius 
turns  shown  in  Figures  4-20,  -22,  and  -24  are  shown  in  Table  4-6. 


Mode 

Pilot 

Velocity 

Radius 

3 

1 

17o7  fps 

300  ft 

1 

2 

17 0 4  fps 

280  ft 

1 

3 

I7o0  fps 

230  ft 

TABLE  4-6„  Course  lb.  Low -Speed  Turns 

From  Section  4„2„4  the  minimum  radius  of  a  coordinated  turn' 
at  17  fps  is  299  feet.  It  is  seen  that  all  pilots  made  tighter 
turns  than  this  (at  17o7  fps  R  325  feet).  However,  none  of 
the  pilots  approached  the  130-foot  radius  recorded  at  18.2  fps 
in  a  57-degree  sideslipping  turn  (Table  4-4). 

The  second  observation  is  that  Pilots  2  and  3  tended  to 
use  roll  control  much  more  than  Pilot  1  in  holding  a  course. 
This  can  be  seen  by  comparing  roll  angle  traces  for  the  three 
pilots. 
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Course  2b  Is  a  high-speed  turn  around  a  promontory  and 
ending  In  a  lOOO-foot-square  area.  The  flying  times  for  the 
various  pilots  and  modes  are  given  In  Table  4-7° 


Pilot  1 

Pilot  2 

Pilot  3 

Mode  1 

220  Sec, 

227  Sec, 

212  Sec, 

Mode  2 

227 

266 

272 

Mode  3 

260 

258 

300* 

Mode  4 

250 

_ - 

210 

254 

TABLE  4-  7.  COURSE  2b,  FLYING  TIMES 


Mainly  this  table  shows  the  added  cautlon--hence  slower  speed-- 
used  with  an  unfamiliar  control  mode.  The  fact  that  Mode  4 
times  were  all  less  than  Mode  3  times  Indicates  the  Increasfe 
of  confidence  with  experience.  The  course  length  for  these 
runs  Is  about  12,500  feet.  Therefore,  a  run  time  of  250 
seconds  Indicates  an  average  velocity  of  about  50  feet  per 
second. 


Figures  4-26;  -2J,  and-  -28  show  the  test  runs  by  Pilots 
1,  2,  and  respectively,  using  Control  Modes  1  and  2,  The 
Mode  1  runs  which  have  headings  dubbed  In  show  and  Interesting 
contrast  In  flying  technique  between  Pilots  1  and  2  and  Pilot  3. 
The  relatively  small  sideslip  angles  of  Pilots  1  and  2  on  the 
straightaway  and  In  slow  turns  Indicate  the  tendency  to  fly  by 
Instruments.  This  type  of  run  requires  fairly  close  attention 


*01d  not  enter  the  target  area. 
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to  keeping  yaw  rate  and  side  velocity  small.  The  smooth  curves 
are  the  result  of  uniform  yaw  rate  and  continuous  side  force 
(due  to  roll).  The  rather  large  sideslip  angles  produced  by 
Pilot  3  show  a  tendency  to  fly  by  '*feel"--that  is,  to  watch 
the  course  on  the  X-Y  plotter  and  control  drastically  when 
a  change  is  required.  The  relatively  slow  vehicle  response 
results  in  rather  large  sideslip  angles  before  a  course  change  is 
is  noticed. 

It  should  be  noted  that  all  three  pilots  were  able  to 
produce  test  runs  as  good  as  those  shown  in  Modes  1  and  2 
at  the  time  these  runs  were  made. 

Figure  4-29  shows  two  Mode  4  runs  by  Pilot  1,  The  first 
run  shows  what  can  happen  to  an  instrument  pilot  when  an 
instrument  is  misread.  The  yaw  angle  meter  has  two  scales --one 
reading  from  minus  180°  to  plus  180°,  the  other  from  minus 
360°  to  plus  360°,  This  is  useful  in  that  it  premits  full 
360-degree  turns  without  saturating  the  meter  and  yet  allows  the 
pilot  to  read  smaller  angles  accurately  by  a  flip  of  the  switch. 
In  Run  1,  the  switch  was  on  the  360-degree  scale  and  the  pilot 
was  reading  the  l80-degree  scale.  When  he  thought  he  was 
approaching  a  90-degree  heading  as  he  rounded  the  promontory, 
he  actually  was  approaching  a  l80-degree  heading.  The  result 
is  obvious.  Run  2  was  made  immediately  to  show  (the  witnesses) 
that  he  was  able  to  control  yaw  as  well  with  the  yaw  pedals 
as  in  the  Mode  1  run  with  the  wheel. 

Figure  4-30  shows  the  "fly  by  feel"  pilot  first  in  an 
aborted  run  which  resulted  from  instinctively  doing  the  wrong 
thing  with  an  unfamiliar  control  and  their  successfully  bringing 
his  craft  into  the  target  area. 
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4 „ 3 o 3 o 4  Obstacle  Courses 

Courses  Ic  and  2c  were  designed  to  test  the  pilot's 
ability  to  fly  to  a  prescribed  point  through  narrow  passages 
and  in  the  presence  of  obstacles „  Figures  4=31^ -32,  and  -33 
show  the  Mode  1  and  3  runs  by  Pilots  1,  2^  and  3,  respectively. 

As  seen  in  the  traces,  the  course  offered  considerable  flexi¬ 
bility  in  choosing  a  path  to  the  dock»  It  is  apparent  also 
that  each  pilot  had  his  own  idea  as  to  the  best  approach. 

However,  only  Pilot  2  was  able  to  consistently  fly  this 
course  without  a  mishap,,  It  is  svaggested  that  in  future 
test  runs  Pilot  2  be  appointed  harbor  pilot  so  that  he  can 
specialize  in  getting  the  open-sea  men  through  the  lOO-foot 
gaps  in  the  breakwater.  One  general  observation  may  be 
made  concerning  low-speed  maneuvering  and  that  is  that  the 
pilots  generally  were  traveling  too  fast  for  the  congested 
harbor  conditions.  At  speeds  below  10  feet  per  second  most 
pilots  were  able  to  negotiate  the  turns  or  at  least  stop 
in  time  to  avert  disasters.  An  estimate  of  speed  is  readily 
obtained  by  measuring  the  distance  between  10-second  blips. 

In  this  way  it  is  seen  that  in  his  first  run  Pilot  1  traversed 
140  to  150  feet  in  10  seconds  as  he  passed  through  the  breakwater. 
Pilot  2  consistently  went  I5  to  20  feet  per  second  outside  the 
harbor  and  10  to  11  feet  per  second  inside  the  harbor.  Pilot  3 
achieved  his  best  runs  at  speeds  of  8  to  10  feet  per  second. 

Another  general  observation  which  is  probably  very  realistic 
for  real  OEM’s  is  that  a  harbor  is  no  place  for  daydreaming  on 
the  part  of  a  pilot.  Many  preliminary  practice  runs  by  the 
three  pilots  included  and  others  ended  much  worse  than  the  first 
and  third  runs  shown  for  Pilot  1,  It  was  generally  conceded  that 
only  by  complete  and  continuous  attention  to  the  controls,  dials, 
and  flight  path  could  a  pilot  successfully  negotiate  this  course. 

t 
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Figures  4-34  through  4-39  show  high-speed  test  runs  in  each 
of  the  four  inodes  by  each  of  the  three  pilots »  The  pilots  were 
instructed  to  fly  at  speeds  in  excess  of  4o  feet  per  second  in 
all  high-speed  runs  except,  of  course,  when  starting  or  stopping. 
Despite  the  higher  speed  this  was  a  considerably  easier  course 
than  the  low -speed  course.  The  strait  has  a  width  of  500  feet 
and  is  the  only  obstacle.  The  open  sea  ahead  of  the  strait 
afforded  ample  room  for  maneuvering  into  a  proper  approach 
course  provided  that  the  pilot  did  not  attain  too  much  speed 
before  starting  his  approach.  The  tendency  of  novice  pilots 
on  this  course  was  to  accelerate  first  and  begin  maneuvering 
afterward.  This  often  resulted  in  speeds  too  high  for 
successful  entry  to  the  harbor  area.  It  is  seen,  however, 
that  all  three  pilots  were  able  to  get  through  the  strait 
reasonably  well  in  all  four  control  modes. 

It  is  seen  in  Figures  4-34  and  4-35  that  Pilot  1  consistently 
made  this  run  at  60  to  70  feet  per  second.  Figures  4-36  and 
4-37  show  that  Pilot  2  made  his  best  runs  at  60  to  70  feet 
per  second  also.  However,  in  Mode  3  he  passed  the  strait 
going  loo  feet  per  second  and  in  Mode  4  achieved  a  similar 
speed  before  braking  drastically  in  order  to  make  his  turn. 

Figure  4-38  shows  that  Pilot  3  had  some  difficulty  in  negotiating 
the  strait  at  speeds  over  60  feet  per  second.  His  more 
cautious  turns  at  20  to  4o  feet  per  second  were  more  successful 
but  were  in  violation  of  the  minimum  speed  Instruction, 
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FIGURE  4-36  COURSE  2c  TEST  RUNS 
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5.0  CONTROL  SYSTEM  DESIGN  CONSIDERATIONS 

5.1  Introduction 

The  main  objective  in  designing  the  control  system  of  the 
GEM  vehicle  is  to  achieve  optimum  performance  of  the  vehlcle„ 

Also  of  great  importance  are  reliability  considerations  which 
often  lead  to  giving  preference  to  a  simple  rather  than  an 
elaborate  system,,  Reliability  can  be  given  more  consideration 
if  control  performance  is  not  critical. 

The  analysis  of  the  inherent  GEM  response  characteristics 
indicate  a  rather  sluggish  behavior  with  respect  to  heave,  pitch, 
roll,  and  yaw  modes.  For  this  reason  the  requirements  on  the 
corresponding  control  components  are  insofar  as  frequency  response 
is  concerned,  easily  satisfied.  In  fact,  it  can  be  shown  that  no 
combination  of  fast-responding  devices  like  electric  transducers, 
transmissions,  and  actuators  with  hydraulic  assist  are  needed 
to  comply  with  the  GEM  performance  characteristics, 

5.2  General  Concept  of  Control  Arrangement 

The  present  GEM  design  concept  considers  two  independent 
power  plants,  one  on  either  side  of  the  vehicle.  The  power 
plants  are  gas  turbines,  and  either  one  drives  via  a  gear  train, 
two  fans  and  one  thrust  propeller.  The  fans  deliver  the  air  flow 
to  the  jets,  while  the  propeller  provides  the  thrust  forces 
necessary  for  forward,  back,  and  yaw  motions.  The  four  fans 
supply  the  air  for  each  of  the  four  base  compartments  and  thus 
enable  the  vehicle  to  pitch  and  roll  by  changing  the  flow  to 
the  respective  compartments.  In  order  to  prevent  severe 
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interaction  of  the  controls,  a  modification  of  the  turbine 
speed  cannot.be  basically  utilized  since  the  power  level  of 
the  fans  and  propellers  would  be  changed  simultaneously.  The 
thrust  control  is  therefore  brought  about  mainly  by  a  pitch 
modification  of  the  propellers  (slanting  the  vehicle  by  pitch 
or  roll  actuation  will  also  contribute  to  the  propulsion  of 
the  GEM  vehicle) . 

The  flow  change  to  the  compartments  is  accomplished  by 
adjusting  the  inlet  guide  vanes  (IGV)  to  the  respective  fans. 

The  hover  height  may  be  selected  by  resetting  the  operation 
point  of  the  centrifugal-force/speed  sensor  of  the  turbines. 
Since  a  change  of  hover  height  is  not  required  very  often,  it 
is  felt  that  the  disturbance  to  the  other  controls  may  be 
tolerated.  '  .  ‘ 

5,3  Control  of  Roll,  Pitch,  and  Heave  Motion 

As  outlined  above,  these  control  functions  are  brought 
about  by  setting  the  IGV  to  obtain  the  desired  flow  to  the 
respective  base  compartment.  The  actuation  of  the  IGV  requires 
positive  positioning  and  represents  (due  to  a  prevailing  high 
level  of  static  friction)  a  severe  control  problem.  Any 
modulating  actuator  would  therefore  require  rather  high  input 
signals  to  develop  sufficient  breaking  force  to  get  control 
motion  starter.  Excessive  overshoot  and  undue  position  error 
may  be  the  result.  To  cope  with  this  situation,  it  is  intended 
to  utilize  either  a  step  actuation  system  (Proposal  AP-5059-R*) 
or  an  air  motor  with  a  self- locking  screw  drive.  To  reduce  the 
effects  of  high  static  friction,  a  surplus  of  actuation  force 
would  have  to  be  provided. 

*Ref .  4 
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In  this  application  where,  in  view  of  the  sluggish 
behavior  of  the  vehicle,  response  is "not  a  major  problem,  it 
seems  possible  to  select  the  gear  ratio  in  such  a  manner  as 
to  provide  (1)  sufficient  force  surplus  (to  overcome  the  effects 
of  static  friction),  (2)  self-locking  action,  and  (3)  sufficient 
frequency  response  with  a  motor  of  relatively  small  size* 


5.H  Control  of  Fore,  Beck,  and  Yaw  Mo^on 

All  of  these  control  functions  are  executed  according  to 
the  present  design  concept  by  modifying  the  propeller  pitch 
angle.  To  attain  fore  end  back  motion  as  well  as  braking 
action,  the  pitch  of  both  propellers  located  in  the  rear  of 
the  vehicle  will  be  altered  simultaneously. 

The  yaw  control  serves  to  change  the  heading  of  the  vehicle. 
For  this  purpose  a  differential  thrust  condition  is  established 
by  increasing  the  propeller  pitch  angle  on  one  side  and  decreasing 
the  pitch  angle  of  the  propeller  located  on  the  other  side  of  the 
vehicle.  At  steady  state,  the  yaw  moment  is  counter  balanced  by 
the  drag  moment.  Because  this  moment  is  a  function  of  speed, 
the  response  of  the  vehicle  to  yaw  input  is  also  speed -dependent. 
At  zero  speed  the  drag  moment  diminishes  and  the  yaw.  control 
approaches  an  unstable  condition.  The  vehicle  has  a  tendency 
to  drift  out  of  heading  position. 

5.5  Description  of  the  Proposed  Pneumatic  Control  System 

Figure  5*1  shows  a  schematic  diagram  of  the  entire  system. 
Compressor  bleed  air  being  used  for  signal  modulation  and 
transmission  is  first  filtered  and  regulated  in  a  common  section 
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for  all  branches  (pressure  regulation  is,  however,  not  required 
to  assure  control  accuracy).  Each  control  branch  represents 
a  pneumatic  bridge  circuit  (see  sketch)  consisting  of  two  fixed 
orifices  (1  and  2),  one  input  modified  orifice  (3),  and  one 
orifice  (4)  adjusted  by  the  INPUT- 

position  of  the  actuator.  Steady-  CONTROLLED 

state  balance  is  attained  when 


FIXED  ORIFICE 


ORIFICE 


Ai/Aa  =  A3/A4.  Thus,  since 
Az/Ai  =  Const,  A4  and,  therefore, 
the  actuator  position  is  pro¬ 
portional  to  the  input,  repre¬ 
sented  by  orifice  area  As.  The 
error  sensor  may  be  a  zero  rate 
diaphragm  and  the  actuator  a 
piston  or  air  motor.  Note  that 
the  signal  transmission  line  in 


(1) 


(3) 


AIR 


SUPPLY 


ACTUATOR 


(2) 


(4) 


ACTUATOR- 

POSITION- 

CONTROLLED 

ORIFICE 


this  all-pneumatic  system  is  tubing 

just  big  enough  to  supply  orifices  of  very  small  size  (approxi¬ 
mately  0,010  inch  diameter). 


A  rather  interesting  component  of  this  all-pneumatic  control 
system  is  the  automatic  trim  control  as  shown  in  Figure  5.2.  The 
error. sensor  is  here- a  pendulum,  free  to  rotate  in  the  X  and  Y 
plane,  A  disc  is  attached  to  the  pendulum  modifying  the  flow 
area  of  four  orifices.  Each  orifice  affects  the  equilibrium 
condition  of  the  control  to  the  IGV  of  one  of  the  four  fans  in 
such  a  manner  as  to  correct  the  flow  to  the  respective  compart¬ 
ments  until  horizontal  orientation  of  the  vehicle  platform  is 
attained.  To  establish  a  horizontal  orientation  initially,  a 
level  Instrument  is  required.  This  instrument  will  be  placed  on 
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FIGURE  5-1  CONTROL  SCHEMATIC 
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FIGURE  5-2 

SEMI-AUTOMATIC  TRIM  PENDULUM 
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the  platform  of  the  hovering  vehicle.  The  hand  controls  for 
pitch  and  roll  will  be  brought  into  neutral  position,  and  then 
the  three  set  screws  of  the  automatic  trim  control  will  be 
adjusted  until  the  reading  of  the  level  instrument  indicate# 
perfect  horizontal  orientation  of  the  GEM  platform.  After 
this  initial  calibration  the  trim  control  will  always  readjust 
the  vehicle  to  horizontal  orientation  whenever  the  automatic 
trim  control  is  actuated.  During  flight  the  trim  control 
will  be  inactive.  In  case  a  new  trim  adjustment  is  required, 
the  vehicle  must  stop  or  drift  with  pitch  and  roll  at  neutral. 
A  short  duration  with  automatic  trim  on  will  then  readjust 
the  level  orientation. 

In  the  inactive  condition,  the  pendulum  is  held  in  trim 
position  by  a  spring-loaded  braking  device.  Actuation  of  the 
automatic  trim  is  accomplished  by  a  pneumatically  controlled  . 
brake  release  (see  Figure  3>2). 
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6,0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  various  areas  of  interest  involved  in  the  studies 
reported  here  necessitate  classifying  the  principal  results 
and  conclusions.  Section  6,1  deals  with  the  simulation  itself. 
This  is  further  divided  into  subsections  dealing  with  actual 
vehicle  characteristics,  the  control  station  simulation,  and 
the  pilot  display  simulation.  Section  6,3  deals  with  the 
control  modes.  Section  6,4  deals  with  the  test  runs  that  were 
used  to  evaluate  the  vehicle,  the  control  mode,  and  the  pilot. 

6. 1  The  Simulation 

6.1.1  Vehicle  Dynamics 

The  simplified  6-degree-of-freedom  vehicle  simulation 
used  here  ha^  been  evolved  over  a  period  of  2  years.  It  has 
been  shown  in  past  reports  that  the  fundamental  analysis  used 
successfully  predicts  the  d5mamic  characteristics  of  both  models 
and  full-scale  vehicles.  The  test  results  in  Section  4,1  further 
verify  that  pitch,  heave,  and  roll  motions  conform  to  those 
found  in  model  tests.  Beyond  this,  the  representations  of  fore 
and  aft  motions,  side  motions,  and  yaw  motions  give  reasonable 
results.  It  is  concluded  that  the  formulation  used  here-- 
includlng  coupling  between  the  various  motions --reliably 
reproduces  the  response  of  these  three  vehicle  motions  to 
external  forces.  This  means  that  the  program  can  now  be  used 
to  evaluate  any  peripheral  jet  vehicle  for  which  the  basic  design 
is  known. 
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The  representations  of  external  forces  are  less  well 
established  than  are  the  equations  of  motion„  The  effect  of 
surface  waves  was  established  in  an  earlier  progress  report. 

However,  waves  were  not  simulated  in  this  study  because  of 

‘  ■  1 

limited  computer  capacity.  The  response  of  the  water  surface 
and  formation  of  a  wake  is  very  reasonably  represented  for 
forward  motions.  However,  it  is  felt  that  a  more  sophisticated 
three-dimensional  representation  can,  be  made  which  will  show  the 
effect  of  rolling  as  well  as  pitching  due  to  the  standing  wave. 

Such  other  surface  conditions  as  a  sudden  change  in  surface 
height,  a  sloping  surface,  or  a  transition  from  overland  to 
over -water  operation  have  not  heen  simulated--although  it  would 
not  be  difficult  to  include  any  or  all  of  these  to  show  how  they 
affect  vehicle  motions. 

The  external  aerodynamic  forces  are  reasonably  represented. 

The  lift,  drag,  and  moment  coefficients  were  based  upon  model 
test  results  where  possible  and  were  otherwise  obtained  analytically. 
The  conclusion  here  is  that  it  would  be  well  to  have  aerodynamic 
model  data  to  verify  the  numbers  obtained  by  analysis.  Two 
parameters  that  are  not  well  known  are  the  moment  coefficient 
in  sideslip  and  the  Interaction  of  side  force  and  drag  during  side¬ 
slip. 


The  characteristics  of  the  air-moving  equipment  (engines, 
fans,  ducts,  and  propellers)  are  well  established.  The  only 
significant  transient  in  this  system  is  the  engine  response  to 
speed-change  commands.  This  command  is  seldom  used  in  practice. 
The  steady-state  gains  can  be  easily  obtained  from  either 
experimental  or  analytical  performance  maps.  The  forces  applied 
to  the  vehicle  due  to  fan-flow  and  propeller-flow  changes  are 
well  represented  in  the  simulation. 
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The  control  components,  which  are  intermediate  between 
the  pilot's  control  and  the  force-producing  components,  may  be 
designed  to  have  whatever  gains  are  required  for  control  of  the 
vehicle.  During  the  test  runs  a  gain  change  was  made  in  any 
control  that  appeared  too  sensitive  or  too  unresponsive.  Initially, 
the  yaw  control  yielded  large  yaw  rates  for  very  small  pilot 
inputs.  The  initially  selected  gain  was  reduced  by  a  factor  of 
5  before  the  Mode  1  runs  were  made.  After  several  trial  runs  in 
Mode  3,  the  gain  of  yaw  moment  to  pedal  motion  was  further 
reduced.  The  conclusion  is  that  the  control  simulation  was 
reasonable  and  that  the  gains  arrived  at  through  studies  similar 
to  this  may  be  specified  for  the  design  of  control  components. 

To  summarize  the  conclusion  of  this  subsection: 

(1)  The  six  equations  of  vehicle  motion  are  valid 
and  can  apply  to  any  peripheral  jet  vehicle, 

(2)  The  surface  effects  of  normal  waves  and  two- 
dimensional  wake  are  well  represented,  but  a 
three-dimensional  wake  representation  and 
possible  variations  from  a  flat  surface  could 
improve  the  surface  simulation, 

(3)  Wind  tunnel  data  may  improve  the  simulation  of 
external  aerodynamics.  However,  the  simulation 
is  quite  reasonable  in  this  respect. 
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(4)  The  air-moving  equipment  is  well  simulated. 
Evaluation  of  any  particular  vehicle  will  require 
performance  maps  and  engine  dynamics, 

(5)  Since  the  control  system  can  be  designed  to 
perform  as  desired,  the  simulation  should  be 
used  to  specify  rather  than  evaluate  the  control 
gains, 

6,1,2  Control  Station 

The  control  station  used  was  very  successful  in  giving 
the  pilot  control  of  basic  vehicle  motions.  The  ease  with  which 
various  modes  of  control--e,g,,  the  manner  in  which  the  pilot 
commands  a  particular  vehicle  motion- -could  be  simulated  is  an 
outstanding  feature.  In  general,  the  human  engineering  of  the 
particular  hardware  used  is  good.  However,  several  limitations 
were  noted  by  the  pilots  during  test  runs: 

(1)  Prolonged  application  of  a  control  that  is 
spring-centered  is  very  tiring.  This  applied 
to  the  push-pull  wheel  motion  and  to  the  wheel 
base  rotation, 

(2)  Controls  that  are  not  spring  centered  should 
have  a  detent  or  other  indication  of  the  neutral 
position.  This  applied  particularly  to  the  yaw 
and  roll  controls. 
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The  pilot  read-out  was  probably  the  weakest  part  of  the 
simulation  in  duplicating  actual  flying  conditions.  Aside 
from  the  absence  of  engine  noise,  salt  spray,  and  wind  noises, 
the  pilot  has  no  feeling  of  motion  and  could  not  see  where  he 
was  going.  Several  Improvements  in  the  pilot  read-out  have 
been  proposed; 

(1)  A  graphic  display  of  heading  would  help 
considerably.  Several  possibilities  have  been 
suggested.  The  simplest  would  be  to  have  a 
small  model  that  would  be  continuously  rotated 
to  the  proper  heading  as  signaled  by  the  com¬ 
puter.  This  would  immediately  tell  the  pilot 
when  and  by  how  much  he  is  sideslipping.  This 
has  a  limitation,  which  occurred  when  the  vehicle 
heading  reached  l8o°--that  the  pilot  must  reverse 
his  thinking  and  yaw  left  to  turn  right.  A  more 
sophisticated  suggestion  is  to  rotate  the  X-Y 
plotter  so  that  the  pilot  feels  his  relative 
rotation  with  respect  to  the  ground, 

(2)  Mounting  the  entire  control  station  on  a  moving 
platform  that  actually  produces  pitch,  heave, 
roll,  and  yaw  motions  will  give  the  pilot  the 
feeling  of  motion  and  response  which  convey  the 
size  of  the  vehicle  to  his  senses. 

Aside  from  these  limitations  and  recommendations,  it  should 
be  pointed  out  that  with  the  pilot  controls  and  read-out  used, 
a  man  can  learn  to  fly  any  specified  course  in  an  hour  or  two  of 
flying  time.  In  shor^  with  the  equipment  presently  in  use,  the 
reference  vehicle  evaluated  is  controllable. 
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6.2  The  Reference  GEM 

The  vehicle  evaluated  during  these  test  runs  was  found 
to  be  stable  in  pitch  heave,  roll,  and  yaw  and  to  be  controllable 
by  means  of  the  fans  and  inlet  guide  vanes  provided.  The  vehicle 
tended  to  be  sluggish,  had  a  tendency  to  sideslip--particularly 
in  cross  winds--and  is  difficult  to  stop  or  turn  to  avoid 
collisions.  However,  the  sluggishness  is  probably  not  unreason¬ 
able  for  a  vehicle  of  this  size  (4o  tons,  60  feet  long).  Also, 
it  was  shown  that  it  can  be  controlled  even  in  a  rather  stiff 
breeze  (40  feet  per  second).  The  pilots  tended  to  feel  a  need 
for  added  side  force,  but  none  of  them  used  the  full  sideslipping 
capability  in  turning.  Yaw  stability  could  be  increased.  However, 
this  is  one  of  the  least  certain  of  the  aerodynamic  parameters 
used  in  the  simulation. 

In  general,  it  can  be  concluded  that  a  vehicle  with  the 
characteristics  of  the  reference  design  can  be  safely  controlled 
by  an  experienced  pilot. 

6.3  Control  Modes 


It  was  found  in  the  test  runs  that  a  pilot  can  learn  to 
fly  the  GEM  with  just  about  any  reasonable  control  mode.  However, 
some  modes  of  control  were  found  easier  than  others.  Despite  the 
early  experience  in  Mode  1,  the  pilots  agreed  that  Mode  4  was 
the  most  pleasant  control  mode.  In  this  mode,  each  of  the  four 
functions  were  separated,  and  no  motion  against  a  spring  was  used. 
In  his  hands  the  pilot  had  only  roll  control.  Propulsion  and 
pitch  were  controlled  by  friction-held  hand  levers,  and  yaw  was 
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controlled  by  differential  foot  pedals.  In  Mode  1,  where  the 
pilot  had  pitchj  roll,  and  yaw  in  his  hands  continuously  and 
was  required  to  hold  two  of  them  against  springs,  the  pilot  felt 
busy  all  the  time.  In  Mode  2,  which  was  the  same  as  Mode  1 
except  that  the  pilot  also  had  continuous  thrust  control  in  his 
foot  pedals,  the  pilot  had  trouble  holding  speed  while  maneuvering. 
He  had  a  strong  tendency  to  accelerate  or  brake  inadvertently 
while  yawing  or  rolling. 

The  recommendations  that  could  be  made  at  this  point  are 
that  the  various  control  Inputs  be  separated  as  much  as  possible 
and  only  those  really  needed  continuously  be  in  the  pilot's  hands, 

6.4  Test  Runs 


In  general,  it  is  felt  that  the  test  runs  used  give  a  good 
indication  of  vehicle  controllability.  It  is  suggested,  however, 
that  certain  other  flying  conditions  be  used  for  pilot  training 
and  perhaps  for  evaluating  the  forces  that  will  act  upon  a  vehicle 
being  designed.  These  added  flying  conditions  are: 

(1)  Docking  in  a  wind, 

(2)  Encountering  a  side  gust, 

(3)  Beaching  across  waves  (breakers) 
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APPENDIX  A 

LINEARIZATION  AND  SIMPLIFICATION 
OF  BASIC  EQUATIONS 

1.0  In  order  to  simplify  future  analysis  of  the  NONR  Reference 

Design  GEM,  the  equations  for  cushion  pressure,  jet  flow  and 

body  motions  as  contained  in  Section  2.4  of  AP-5047-R,  will 

be  linearized  and  written  in  nondimens ional  forms . 


1.1  The  general  nonlinear  transient  pressure  equation  (Equa¬ 
tion  106,  AP-5047-R)  is: 

•  P  p 

Pi  “  •  dr  (Qi  -  2_  ‘^iK^  -  (-  *  +  -  Vi})  (A-i) 

^  lUl 

It  is  assumed  that  if: 

1.  The  base  area  (Ag)  is  divided  into  4  equal  areas 
(A^);  Aj^  =  Ag/4  • 

2.  The  projected  jet  length  (b)  in  the  fore -aft  direc¬ 
tion  (x)  is  twice  the  length  (a)  in  the  port- 
starboard  direction  (y);  b/a  »  2 

3*  The  quantity  ^  as  defined  in  Section  2.4  (AP-5047-R) 
is  equal  to  0.5;  =  0.5 

Then  it  can  then  be  seen  that  ' 


_  ba  b 

^  -  K^'  2  (a  +  b)  -  3 

and  a  time  constant  xa  can  then  be  defined  aa: 


(A-2) 
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Dividing  the  pressure  equation  (Equation  A-1)  abov*  by 
W^/Ag  and  multiplying  it  by  yields: 

bV.  A  . 

3g-Jx  (57^)  Pi 

.  Ap  •  _  •  _  •  /*\ 


bV.  A„  P„.  ,  A 

3FRT  RX7  ^^i  ■  Z 

^  ^  K»1 


From  the  perfect  gas  law^  the  total  pressure  in  a  com¬ 
partment  may  be  defined  as: 


where  p  has  dimensions 

‘’ll  '  P  *  **  ..  lb-8.C.* 

— rn — 


substituting  above  and  A^^  =  Ag/4  into  Equation  A-3  yields: 


n 

Z  Oik) 


lUl 


bV i  P^n  •  _  •  • 

- 3  (  -  z  +  (A-4) 


By  dividing  by  and  by  in  Equation  A-4 

it  can  now  be  written  as: 


CiG^ 


/  ^i  .  QfiK  QjKv 
^  ^fiK 

y^i 

n - tr^ 
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Considering  the  effect  of  a  wave  moving  beneath  the 
vehicle  as  In  Equation  (107)»  AP-5047-R. 


n 


Ti  (gr)  =  ■ 


3^  "j'  ^ 


flK  '^IK 


flK 


) 


bpV^A 


+  if  (^)  cos  ^  (vt  +  x^)j 


since  Kj^  =  0.5 


b 

5 


C 


therefore 


4Cj^b 


^  =  g  (a  t  =  2  A. 

3  ^  B 


G 


above  pressure  equation  becomes 

A  ^  ^ 


'  V  "I'-"'  -W 

^T^B 


W. 


z  *1*  yi* 
E  +  TT®  ■  TT^ 


w  X  +  *i) 
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Appendix  A  of  AP-50^7-R  contains  a  simplified  form 
flow  equations. 

Equation  {2C))  in  Appendix  A  (AP-5047-*-)  !• 


Qi 

^fi 


with  F 


^GT  ■  I  ®i 


MS  -  5  sin  6^ ) 


-  ^i  ®: 


j 


'‘G  ■  5 


®lo  “  ^ 


^  ■  5  pVp  ■  ‘  5 


Q 


Multiplying  Equation  (A-l2)  by  tt; 


Qi  2pV/^BG 


R  ^Gj  '  I  ^  pvp 


-  (g  -  I  sin  ej)  py^j 
?3|[  “  [^TT  ■  5  E  ^  Pi 


-  (1  -  I  g  sin  0.)  g5 


B 


') 


Pit 

where  hj^  *  h  -  z  +  Xj^0  r^d  +  sin  ^  (vt  +  Xj^) 


h  being  the  steady-state  operating  height 


of  the 

(A-9) 

(A-10) 

(A.11) 

(A-12) 

(A^13) 

(A.U) 

(Ar-15) 
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Equation  (iUl9)  is  applicable  for  a  range  given  by 


-PV.^G 


^  (P,  -  P,)  , 


For  the  third  condition  on  the  compartment  jet  the  follow¬ 
ing  relationship  holds: 

^Ki  "  ^fiK  ‘  ^iK 

using  the  value  given  in  Equation  (A-16)  for  an  expression 
for  is  obtained: 


‘iK 


fiK 


=  Fa 


1  -  Fi 


‘iK  ^^i  "  ' 


CA-20) 


Equation  (A-20)  is  applicable  for  the  range  as  given  by 


^Pi  ^  ^iK  ^^i  ■ 


GiK  2pV  *A 

Ut  TTs  =  -jr-  TTa  ^ - -inf 

1  < 


apVj'A^ 
C 


Multiplying  Equations  (A-16),  (A-I9),  and  (A-20)  by  ir* 
and  redefining  the  boundary  conditions  yields: 


TTa 


^iK  ^iK 


iK 


7f^  “  TTa  -  Fa 
^i  ^fiK  ^i 


Ci h.,,  AfiPj 


ir  J.  9V  ''iK  "iK  /"B**i  ^B^kJ 

—  TTa  +  ZFi  ft —  c —  {-rj — - D — Jl 

H  ^i  **i  *'c  ''c  J 
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when  j—  (het - g— 

1  c  c 


-5F* 


^iK  ^IK  1  ,  ^iK 


^iK  ^iK  ^B^K> 


’IK  ^IK  I  _  IK.  p  IK.  IK.  t  o  .  o 


(A-22> 


^2  ^  /^B^l  ^B^Kx  .  .  JLa 

when  -  ^  ^  ^  (.^  -  ^ 


_  ^IK  ^iK 


=  Fs 


^iK  ^iK  /^B^l 


IK.  _  OP  IN  IN  /  P  *• 

■p —  -  2^1  p —  yt^  ^nj 

Cl  Cl  Hi  w 


A  P 
**c 


fA-23^ 


TT  If  ^R^lf 

,,u-_  ,  tr2  .  IK  /  B  1  B  K\ 

when  +  sjr-  15^  (-g- - g— ) 


1  -  ,  ^*1K  -n  TT  ,  ^iK  1  TTb  ^^1K\ 

1  -  z  +  -T—  B  -  -— —  (i)  +  -r-  sin  Y" 


(A-2l^) 


Equation  (35),  AP-5047-R,  is  the  general,  nonlinear  ex¬ 
pression  for  heave  orations. 


tt  n  r  **  n2  ’li.i 

T  *  “  ■  Z  '  P  r  Z  CJCJ  ® J  Z  Z 

*  1=1  c  i=l  ^  ^  i*l  K*l-» 


i-1  K.l+1 


^fiK  ]  yf 

^iK^iiy  ‘ 


At  this  point,  if  the  base  pressure,  Pi,  is  thought  of  as 
being  the  average  steady  state  value  plus  a  variation,  it  can 
be  seen  from  the  above  expression  that  at  steady  state  condi- 
tions,  i.e.,  z  =  0,  the  average  steady  state  pressure  is  less 
than  the  quantity  W^/Ag  by  an  amount  that  is  equivalent  to  tha 

thrust  produced  by  the  Jets,  or; 
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nMWl,.<nMW 


r  n 


Xi  = 


^B^AV 

"177“ 


=  1  - 


^  sin  e. 


i  i.l 


IT 


n-l 

1=1  K-i+1 


-i+1 


(A-25) 


Defining  a  dimensionless.,  normalized  pressure,  as; 


^B^i 


-  Xi 


(A-26) 


since  Xx  is  constant. 
Pi 


^B^i 


dimensionless  body  motions; 


*  “  K 


-5 


be 

5E 


and  dimensionless  flow; 

Let 

Let  T2  ■  -»  T  1 


(A-27) 


fA-28) 

(A-29) 

rA-30^ 


(A-31) 


Equation  (A-8)  may  now  be  written  as: 
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Atituu 


T  1  = 


TTi  T3i  +  TTa  Y.  ^  ^iK 
K=1  *■ 


“  T2  -  Z  +  --g- 


^  ,  2*1  .  2y^ 


^  /27rv.  TTb  /2vt  ^ 

+  TT  Xi  (-3-  +  ^ 


'A-32 


Note  here  that  X  in  the  expression  for  wave  form  is  now  Xz 
to  differentiate  It  from  the  pressure  expression  Xi. 

Dividing  Equation  (A- 14)  by  the^ steady-state  fan  flow 
(Qflo)  yields: 


TT:  =  2Kj^[  (tT  -  I  h  Pi  ■  ■  5  H  wf  Plo]  ^fl 


where  = 

and  = 


Oi 

^cu, 


Linearizing  on  the  basis  thst 


Pi  =  Pi,  +  Pi 


(steady  state  plus  a  variation) 


2.  Fan  flow  consists  of  a  steady  state  plus  a  variation. 


From  Equation  (A-I5) 

»  h  -  2  +  Xj^e  -  y^0  +  Ej^  sin  (vt  +  x^) 
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Afimnt 


1 

5 


sin  6. 


'j]sf  (•’l  +  Pl> 


(1  - 


1 

5 


sin 


>4 


bTTx'Q. 


^fic 


BTi 


fl 


^  ‘Ifi 


(A-33) 


From  Equation  (41)  (AP-5047"R.) 

BQ. 

aq^  '5i  -  2Pi 

*51 

therefore  _ 

BTTlQf 

-  ^Pj 

SQfl 

=  ■  [^TT  ■  5  E  Pi  ■  5  E  wf  ^io] 

(A-34) 
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Equation  (A-35)  may  be  written: 

=  2Kj^  ^  (1  -  I  g  sin  ep  (Xi  -  -^y^)  " 


•j(Xi  + 


[■ 


+  Xi  -  2  + 


jlfl 

2x, 

a 

% 

TT 


1  G 
S  K 


sin  6 


'i>] 


Trb  /2vt 
Xi 


2x 

+  -xi) 


] 


+  (1  ■  5  H  ^  ^1 


^  +  *lfi)  ‘Ifi 


[- 


‘i  -H  ^^i  —  ^i  j  Tb  / 

2  +  -^  -g  -  _  e  +  -jj.  sin  ^  ( 


2x^ 

“F" 


2vt 


2x 

-IT 


(A-36) 


It  is  noted  here  that: 

\  *  *Hca  ^  (positive) 

“  \b  when  TT^  15^  ^  0  •  (negative) 

Equations  (A-2l),  (A-22),  and  (A-23)  may  be  rewritten  by 
setting 


^B^i 

nrr 


Pi  +  >-a 
P|^  +  Xi 
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TT  K 

irr  =  *^Lik 


Tlk 


(A-37) 


Boundary  conditions  for  Equation  (A-37) 


hik  -  (1  - 


F2 )  ^"2 


K, 


Pik 


''ik 
SFaF  1  7^ — 

^i 


when  Fi  (p^  '  Pk^  ^  I 


1  ^ik 
^Uk  “  5 


"Pik 


\  X  ^Xlt  ^  “  1 

/  when  -  ^  ^2  ^  Fi  (pj^  -  Pj^)  ^  +  5 


^Ik 

*4.ik  = 


Kjik  -  2F.F, 


1  —  * 
when  2  TTa  ^  Fi  jjj*  (p^  t  pj^) 


TTa 


AP-306I-R 
Appendix  A 
Page  13 


fliResearch  Manufecturing  Dvision 


Pk0tiUM, 


1.2  Equations  (35),  (36),  and  (37)  of  AP-5047-R  present  the 
equations  of  motion. 

Heave 


AiPi 


fi 


sin  d 


j 


•  T  i 

1=1  kFi+1 


Q'fik 


+  W 


c 


Revnrlttlng  the  above  equation  using  the  notation  as  expressed 
In  Equations  (A-25)  and  (A-26) 


(Ap-38) 


where  =  %/g7H 


Roll 


•• 

T  ^ 


n 


Y.  ■ 

1=1 


I-  n 

1-1 


n-1 


n 


sin  + 


I  I 


^Ik^^flk 


1=1  k3.+l 


With  Equations  (A-25)  and  (A-26) the  above  equation  may  be 
\7rltten  as  follows : 


(A-39) 
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Pitch 


f®  =  Z  \  Vi  +  p  [Z  ‘1"  *  Z  Z 

1=1  ‘-1=1  1  1  1^.1  v.^1. 


^ikQjik 

^ik^ik 


1=1  LUl  ^  ^  ■  1=1  k=i+l 

Equations  (A-25)  and  (A-26)  into  above  equation  yields: 


] 


ee 

^  = 


2^1- 

nr  Pi 


(A-40) 


where 


AP-506I-R 
Appendix  A 
Page  15 


^iPeseardi  MamifiH:turin9  Division 


Tl  Ps  *  -  ■JTi  ^3  +  7r2  Q32  +  Qs4) 


r  /2tv  ^  Trb  /2vt  ^AlL9^ 


P4  =  -  TTi  ^^4  +  TTa  ^543  +  ^  IJ*  1 ) 

-T2[-i  -  5‘9  -  -ji  +  H  ‘=°® 


Trb  /2vt  1x1  rA_441 

Xi  ^-F"  ■  J  '  ' 


2 .2  Flow  Equations  f Peripheral  Jets^ 


Applying  Equation  (A-36'  to  each  sector,  neglecting  the 
products  of  the  variations,  p^,  z,  etc.  Reasons  for,  and 
proof  of  validity  of  this  linearization  vd.ll  be  given  In 


Section  2.5. 

TTi  ?5i  5=  '  2  E  ^ 


5  .  Ijs  j  ^  .1„  g  (Svt  .  2|,,j 

•  (1  ■  5  K  ® j )  (^1  +  y 


+  2Kka  (  1  -  5  E  Sin  Sj)  (X: 
TTiJ^a  •  2Kj^^  -^(1  ■  5  ^  8ln  6j  )  pa 


A«P 

i  -  -^) 


+  Xi 


TTb  /2vt 


^)] 


-  (1  -  I  ^  sin  9j)  {X.J  +  -J^) 

+  ^*HcA  ^  ^  ^  E  ® i  ^ 

c 
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Ariam* 


TTi  Qs  =  2Kj^  j^(l  -  5  p  «ln  Sj  )  ps 

^  X.  [-  I  .  2^  ^  ^  ^  Sin  g  (iji  ^  ^)] 

-  (^  ■  "g  ^  si*'  ^ j  )  {^1  +  ^)  ‘Ifjl 

+  ^  ^  ’  5  E  ® i  ^ 

c 

TT:  ^4  =  2Kj^^  {(1  -  ^  g  sin  Sj)  P4 

+  X.[-i  (2vt  .  gxajj 

*  (1  *  2  p  sin  @j)  {Xi  +  °)  q£^j 

+  2Kj^  ’  5  E  ^ 

c 


\i  =  •SCA  if  TTi  75^  is  positive 
=  Kj^g  is  TTi  is  negative 
2.3  Compartment  Jets 


Applying  Equation  (A-37)  to  each  compartment  Jet,  re 
taining  the  products  of  variations,  pj^,  s,  etc,: 
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C12 

rr 

T^a  15-2 

■  ''l.=  - 

tpjj  (pi  - 

Pz)  Hia  ^ 

C  12 

irr 

TTz  Z52I 

-  ^-2  TT 

=  Cl- 

-§r--'5 

H 12 

=  1  -  Z 

+  J  + 

sin  I| 

2vt  1 

T* 

C23 

Cl" 

^2  ^23 

- 

"P.,  ■ 

Ps)  1^23 

C23 

cr 

^2  T332 

'§f 

■  §f^  j 

H23 

=  1-2 

+  ^-3  + 

J 

^2  ^34  =  Kl3^  -  (P3  -  P4)  H34  ^ 

^2  ^43  =  ~  ^43  ^ 

TT  1  3  ^  A  j.  ^34  -4—  2vt 

h34  *  1  -  z  -  ^  0  +  Sin  ^  -g-  J 


^  TTg  154:  =  -  Kp^^  (P4  -  pa)  H41 

TTs  ^14  =  ^  Wz  -  ^  TTz  ?j4a 

1i4a  .  I  .  I  -  1  ^  Sin  g  (^  -  |), 


XiK 


"PIK 


2  (1  -  Fz)  ^ 
2  Fa  Fi  pi* 


Vz 


"IK 


If  Fa;;^^H 


"IK  ^Pi 


(A-49) 


(A-56) 


(A-5I) 


(A-52) 


-  Pr)  i  -  5 
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K 


LiK 


TiK 


XiK 


1  ^  ^IK 
Fi  ^ 


-  5  "^2  C“  ^  Cp  ^iK  ^Pl  •  Pk^  ^  5 


1  ^IK 
2  Fa  2  TTa  ^ 


K, 


PiK 


2  Fa  Fi  g — 


iK 

i 


^  >if  5  ^2  ^  ^  Fi  (p^  -  pj^) 


2.4  Body  Motions 


Equations  (A-38),  (A-39),  and  (A^40)  become: 


Heave 


—  a  ^  ^  ^ 

*  “  -  tSj  (Pi  +  pa  +  Ps  +  P4) 


tA-53) 


Roll 


W_a 


^  BT —  (‘pi  ■  Pa  +  Ps  +  P4) 


(A-54) 


Pitch 


g  ,  _  _  _  _ 

Ij^f  “BX^  {-  Pi  +  Pa  +  Ps  ”  P4) 


(A-55) 
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2.5  As  previously  noted,  the  products  of  the  variations 
(Pl,  etc.)  were  neglected  In  the  peripheral  Jet  flow 

equations  and  retained  In  the  compartmentatlon  Jet  flow 
equations. 

Consider  the  product  of  two  variables  X  and  Y,  with 
steady-state  values  of  Xo  and  Yo.  If  x  and  y  are  variations 
froQ  steady  state,  then 

X  -  Xo  +  X 

Y  -  Yo  +  y 

and  the  product  XY  »  (X©  +  x) (Y©  +  y) 

-  XqYo  +  Xoy  +  YoX  +  xy 

Let  X  here  denote  a  dlrnenslonless  pressure  term  and  Y  anQther 
dimensionless  variable. 

2.5.1  For  the  peripheral  Jet  Xo  •  1,  x  «  0.1 

Yo  -  1,  y  -  0.1 

therefore.  Including  the  product  of  the  variation 

X  -  1.1  Xo 

Y  -  1.1  Yo 

and 

XY  -  (1  +  .1  +  .1  +.01)  -  1.21 
neglecting  the  product  of  the  variation 
(1  +  .1  +  .1)  -  1.2 
The  error  Is  seen  to  be  O.83  percent. 
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fkttiUM,  Alumt 

For  the  compartmentatlon  jet,  Xo  -  0,  x  >  .1 

Yo  -  1,  y  -  .1 

Including  the  product  of  the  variation 

XoYo  +  Xoy  +  Yox  +xy=0+0+  .1+  .01  *  .11 

and  neglecting  the  product  of  the  variation 

XoYo  +  Xoy  +  Yox  =0+0+  .1*  .1 

The  error  is  seen  to  be  10  percent. 

Therefore,  the  product  of  the  variations  were  re¬ 
tained  in  the  compartmentatlon  flow  equations  and 
neglected  in  the  peripheral  jet  flow  equations. 
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3.0  Completely  Linear  Equations 


Let: 


-  ^2’ 


(A-56) 


(A-57) 


2Kj^  {1-1/2  ^  Sin  fj)  (Xi+  (A-58) 


Kj,  =  2Kjj^  (1-1/2  g  Sin  Sj) 


(A-59) 


Kg  =2Kg^Xx 


(A-60) 


r  A-P^ 


^  ^io  =  ^hA  E  Sin  ej)(X,  -  (A-61) 
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3.1  Peripheral  Jet  Flow  Equations 

The  above  substitutions  Into  Equations  (A.43),  (A-46), 
(A-47)^  and  (A-48)  yield: 

TTl  IJl  =  TTl  Ql  0  +  Pi  -  Kj) 

-Kg  [z  +  B  e  -  a  3  -  ^  Sin  (^-  B)]J 

TTl  Tja  =  Tl  ^2  0  + 

-Kg  [i  -  B 

TTl  Qa  =  TTl  Q30  + 

-1^  [z  -  B  ?  *  I  ?  -  ^  sin  3$  (2g£  *  B)  j 

Vl  Q4  =  TTl  Q40  + 

-Kg[z  +  B0+a?-^Sln5|  -  H)]J 


(kc  P-  -  %  ^2 

?  -  I  3  -  ^  sin  +  B)]j 

{kc  P3  -  h  ^f3 


.(A-6e) 


CA-63) 


(A.64) 


.(A-65) 
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LET:  5  Tra  ^  »  5  Ta  ^  -  1^ 

"*0 

^  ^ 


In  Equations  (A-41),  (A-42),  (A-43),  and  (A-44),  let  -  0 
and  In  the  perlpherial  jet  flow  equations  let  q^^  > 

Q  4.  ^  D  Q  +  ^  jfl  r 


Substitution  of  the  appropriate  flow  terms  from  equations 
(A -62)  through  (A-69)  into  each  pressure  equation  yields  after 
cancellation  of  steady  state  terms  (iti  a  Kp  +  Kq) 

(Kg  +  -  Kp  -^-)  Pi  -  Ky  Pa  -  P4 

=  ^eC*  +  B?-a0-|sin^  -  ^)]  +  Kp 

+  Ta[*+|t  -  -  |  (^)  cos  ^  ■  5O  (A-7O) 
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fkttnU*  Arim^ 

^  S^p?  ^  -  Ky  Pi  -  Ps 

=  Kg[z  -  B  e-  a  3-  ^  Sin  (^  +  B)]+  Kp  q^g 

+  T2tz  -|  e  -  I  i  ^  (|p)  cos  g  (^  +  |)  ]  (a-71) 

Sgf"  )  Pa  -  Kjj  Pa  -  Ky  P4 

=  Kg[z  -  B  e  +  i  ?-  E  ^  ^  H)]  +  Kp  qf3 

=  T2[  t  -  5  e  +  n  -  g  (^)  Cos  g  (^  +  1)]  (A-72) 

*Si  •  *D  aTT  )  P*  -  *Sc  Pa 

=  Kg  [  2  +  B  e  +  i  ^  -  E  Sin  g  (2^  .  'B)]  + 

=  t^55-|(^)Cos5|(^.  1)]  (a.73) 
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3.3  Solution  of  Pressure  Equations 

3t3.1  Dividing  each  pressure  equation  by  the  coefficient  of 


dn 


^  * 

Pi'  (K  +  K  +  K  ■  K j  - -),  and  denoting  the  right  sides  as 

,  Equations  (A-70),  (A-71),  (A-72),  and  (A-73)  are  now 
written  in  matrix  form. 


'  1 

-B 

0 

-A~ 

- 1 

"Ol 

Cl 

-B 

1 

-A 

0 

P* 

Cs 

0 

•A 

1 

-B 

• 

Pa 

Cs 

-A 

0 

-B 

1 

_ 

_P4 

C4 

(A-74) 


where 


A  » 


K. 


Kc  +  K,  +  Ky  -  Kj 


B 


K. 


Kc  +  ■  *d  - 
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The  general  solution  of  Equation  (A-7^)  !>: 


B(1^A*-B*)  2AB  A(l+B*-A®]  Ci 

,  B(l4A*-B*)  (1-A*-B®)  A(l+B*-A*)  2AB  Ca 

*  _  -  _  -  • 

1-2(A®+B®)+(A®-B*)®  2AB  A(l+B®-A®)  (1-A®-B®)  B(l+A®-B®)  Cs 

A(l+B®-A®)  2AB  B(l+A®-B®)  (1-A®-B®)  C4 


from  Equations  (A-70),  (A-71)i  (A-72),  and  (A-73)> 
can  be  written  in  the  general  form; 

^  ^  ®  ^  ■  "H  ^  ^  ^fi 


(A-75) 


+  T2  [?  ±  -5  i  ±  -5  -  I  {^)  cos  ^  (2^^  {^)±  |)] 

where  the  signs  of  the  “e,  ~4>,  and  b  terms  depend  upon  which 
pressure  equation  is  being  expanded. 


Completing  the  solution  of  the  pressure  equations 
now  proceeds  with  the  substitution  of  the  proper  terna 


in  the 


[c.] 


matrix. 


3.3.2  Flow  Terms 

The  Or.  terms  can  be  written  as 


[c,]- 


Kc  VK,  +  -  Kd 


•N 


(A-76) 
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since  all  of  the  z  terms  are  positive  and  the  matrix  in  Equation 
(A-75)  Is  symmetrical,  the  coefficient  becomes.: 

1  -  A»  -  +  B  +  BA®  -  B®  +  2AB  +  A  +  AB®  - 


or 

(1  +  A  +  B)[l  -  (A  -  B)®j 


AP-5O6I-R 
Appendix  A 
Page  30 


<1iRM«arcli  Manu^Kturing  Division 


Wiwit,  AHtm» 


3.3.4  Pitch  Terms 

The  pitch  terms  in  may  be  written  as  (denoting  ^  as  s7) 


Kg  'B  (1  + 


2Kg^ 


•) 


fl 


4l 

-1 

-1 


(A-79) 


Assigning  the  correct  sign  on  each  coefficient 
yields : 


on  Pi  :  (l.A®-B®)  -B  {1+A®-B*)  -  2AB  +  A(l+B®-A*), 
on  pj,  :-(1-A2.B®)  +B  (l-fA®-B®)  +  2AB  -  A(l+B*-A*), 
on  Pa  :-(l-A®-B®)  +B  (l-fA^-B®)  +  2AB  -  A(l+B*-A*), 
on  p^  :  (1-A*-B*)  -B  (l-tA®-B®)  -  2AB  -  A(l+B®-A*). 

The  P2  and  pa  terms  are  of  opposite  sign  then 
the  Pi  and  p^  terms,  therefore 


1  -  A»  -  B®  -  B  -  BA®  +  B®  -  2AB  A  +  AB®  -  A® 
or  (1  +  A  -  B)  [l  -  CA  +  B)®] 

Ti 


£i 

5. 


l-2(A*+B*)+(A*.B®)® 


^^fl 

dTl 


-1 


-1 


•  ) 


(A-80) 
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3.3.5  Roll  Terms 

Using  similar  notation  as  In  pitch,  the  roll 


term  In 


may  be  written  as: 


Kg  a  ( 1  +  *  _  ) 

2Kg  a 


^fl 

a? 


<t> 


(A-81) 


Assigning  the  correct  sign  to  each  coefficient 
yields : 

on  :  -(1-A®-B®)  -B(l+A®-B®)  +  2AB  +  A(l+B*-A*)^ 

on  pe  :  -(1-A®-B*)  -B(l-fA®-B®)  +  2AB  +  A(l+B*-A*)^ 

on  Pa  :  +(1-A*-B=)  +B(1+A*-B®)  -  2AB  -  A(l+B*-A*)^ 

on  p^  :  +(1-a2-B®)  +B(l4A*-B*)  -  2AB  -  A(l+B*-A*). 

The  p^  and  ps  terms  are  of  opposite  sign  then  the 
Pg  and  P4  terms,  therefore 


1  -  A®  -  B*  +  B  +  BA®  -  B*  -  2AB  -  A  -  AB®  +  A* 
or  (1  -  A  +  B)  [1  -  (A  +  B)®] 


ffl' 

(l-A+B) 

[WA+B)»] 

L  ♦. 

1-2(A«+B®)+(A®-B®)* 

r 


Kc+K^^4Cy-Kd 


^fi 

afi 


1 
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3.3.6  Wave  Motion  Term< 

Equations  (A-70)  and  (A>73)  contain. 


*  •'e  I  fl  ■  '*> 

or 

“  I  ^ 

“  I  ^^Hcos  ^t  cos  +  sin  igt  sln^) 
Similarly,  Equations  (A-71)  and  (A-7S)  contain, 

-■^1 

■  I  5?;  ■  ^  ‘  siS' 
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The  wave  terms  in  C4  I  become: 


M 


['*]■ 


«•  E  _j_  2tv  _ _ Trbb 

~  Ti  ~Ti  ^  ~ 


-  T2 


I  (^)  cos  t  CO. 


jy 


'C  +  *X  +  *Y  -  *D  ^ 


a<i 


fi 


1 

1 

1 

1 


'{A-83) 


+  i 


+  K. 


2ttv 


Trbb  ^  E  f2iTV\ 


E  Ti  °°®  "Ti  ■  ■'2  Ti  (-5Ci) 


g-TTV 


t  sin 


3Qfi 

Kc.Kx.K,.K,^ 


1 

-1 

-1 

1 


Assigning  the  proper  sign  to  the  coefficients  yields: 


t.,  V  E  „  4  _  27rv  _ _ Tfbb  2t'vX’s  E _ Sttv  ^  Trb 

l-(A-B)®]  T*^E  1  "Ti  -T^  - 

-  l-2(A«+B»)T(A*iB^’^ - 


3 

ST, 


fi 


[l-jA.Bl  j  k  I  -i"  ^ 

r-2(A«+B«*^)+iA*-B*)*  * 

If  .  IT  .  IT  tr  rr 


*c  - 


3 


1 

1 

1 

1 

+1 

-1 

-1 

+1 


(A-84) 
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4,0  Equations  of  Motion 

To  complete  the  linearization,  the  expressions  ^r  the 
transient  pressure  in  each  compartment.  Equations  (A-70),  (A-71) 
(A-72),  and  (A-73)  will  be  combined  with  the  equations  of  motion 
Equations  (A-53),  (A-54),  and  (A-55). 

4.1  Heave 

From  Equation  (A-53) 


where  oijj  =v/g7K“ 

Substitution  of  the  pressure  terms,  as  given  in  general 
form  in  Equation  (A-75)>  into  the  above  equation  involves 
inclusion  of  the  terms  defined  as  the  ^  matrix  which 
have  an  effect  on  the  heave  motion.  This  formulation, 
which  deals  with  a  symmetrical  platform  will  therefore 
exclude  pitch  and  roll  terms  in  the  heave  mode  equations. 


AP-5061-R 
Appendix  A 
Page  35 


^il^search  Manufacturing  Division 


The  applicable  terms  are: 

q££  Equation  (A-76) 

s  Equation  (A-78) 

Wave  form  Equation  (A-84) 

•e 

Denoting  z  as  and  making  the  above  substltu 

tions  yields: 


-  (1  +  A  +  B)  [1  -  (A  -  B)®]  Kg 


1  -  2(A*  +  B*)  +  (A®  .  B®)*]  ^  Kc  +  ^ 


(' 


+  KJ  +  51^  (if,  +  Ifj  +  If,  +  if,) 


'I 


/E  2-irv  _ _ Fb5  .  2tvt2E _ 2Trv  _ _  irb  I 

(j.  sin  _  t  cos  -^  +  cos  -5;;^  t  cos 


Note  that: 

1  -  2  (A^+B^)  +  (A^-B®)*  -  {1+A+B)(1+A-B)(1-A+B)(1-A-B) 
and 

(1+A+B)  [1  -  (A-B)2]  =  (1+A+B)[l  -  (A+B)].[  1+  (A+B)] 
Therefore,  the  coefficient  becomes: 


(1  -  A  -B)  (Kc  +  K,  +  Ky  -Kd^ 
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Equation  (A-83)  becomes: 


r  ApP  j*l 

j)  [1  -  +  -If)  -g^] 


(1  -  5  ^  sin  e 


ArP 

(Xi  +  y  °). 
"c 


« [>  •  I- •  >  %1 


*^f2  ‘^fs  ■*” 


Xi  E 

TTS 


(1  -  ^  g  .In  e,)  h 

*  fi  rx  ,Vo,^fil 
[1  -  (xi  +-„^) 


2irv  ^  TTbB 
sin  ^  t  cos  -5^ 


.  2irvT2 


E  2tv  ^  _  ■'Tb 

TT  ' 
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4fimn§ 

(1  +  A-b{i  -  A+B)^] 
1-2(A®+B®)  + 


BQ 


1  -  A+B  = 


•'c  -  % 


21V 


Kc 


"STJ 


fi 


Equation  (A-87)  becomes: 


5 

.W3*/w 


fi 


?  = 


c+V*Sr-KD  wr 


2  ^'<ff 


(1  8)  1?  + 

2KeB 


fi 


*^fi  *^fg  ^fs  ~  ^f4 

4 


Kc  +  *Sc  ^  ‘V  - 


r  _ 2ttv  ^  _j_  TTbS 


cos  ^  t  sin 


-E-  ^ 


.In  Is  t  sin 


TTb  T 


(A-88) 
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AHmut 


Coefficient  of  ^  terns: 


-  V 


2K^X,  H 


2K„.  (1 


-  I  g  Sin  S,)  1  -(Xi  2Fi  §44 


Xi  ^ 


M  1  G  .  vf,  ^B^o.^^fll  a  Fi 

(1  -  ^  ^  sin  -(a:  +  ^ 
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Coefficient  of  q£^  terms: 


^^fl  ~ 


2Kj^  (1  -  5  ^  sin  )  (\i  +  °) 


•  5  E  “§f^j  ^ 


(1  -  ^  g  sin  Sj)  {X-  +  -§^) 


(1  -  ^  g  sin  e^.)  1  -  (xa  +  -^)  -g|i  j  +  ^  ^ 


Coefficient  of  wave  terms: 

K  E 
^  E 

d^JfT 

^  (f)  _ 

*''kA  ■  5  K  •^"  *J>[  1  -  (>•!  +  -^)  -g^]+ 

_ _ ^ 

(1  -  If  aln  «j)  [l  -  (>.i  +  -^)  +  i$B  |^~ 
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Equation  (A-88)  beconas: 


(1  -  1  g  Sin  ®j)  [l  -  (Xi  +  nj^) 

]  ^  '5?b)  ^ 

(1  “  -g  ^  sin.  6j )  (Xi  +  ) 

(1  -  )§  sin  ej)  [l  -  (Xi  +-5^)  -^] 

i  (qf ,  -  qf^  -  qf3  +  qf J 

(1  ■  E  H  ^ 

{Sttv  _  ^bl5  2irvxz  E  giry  ^  ^rb  I 

CO.  -J-  t  .in  ^  -  —5^  x,E  “51  ‘  SXlj 


(A-89) 
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Coefficient  of  $  terms: 


(A-90) 


I 


-  2Kj^  Xa  a 


(1  -  5  H  ■^"  [^  •  ^]  +  <S+S>  1^ 
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Coefficient  of  tenes; 


Equation  (A>90)  becomee: 


^^fi  ^fa  ’  **fa  "  ^£4)  (A-91) 
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5.0  Recapitulation  of  Equations  of  Motion 
5. 1  Heave 

Equation  (A-86)  may  be  written: 


[ 

r 

r 


•*  +  2^2  ^ 


21  z 


^  2  E  2;nv  ^  ^  2TrvT2  _ E_ 

+  tr  aln  — t  cos  — — 

K 


Sttv 

“z"  1  r"-"  —  t;  -5c3r  “>*  -rr 


? '  -  aSl 


(A-92) 


In  Equation  (A-92) 


1 

[ 

[ 


(PZ 


(1  -  J  I  .in  8j)  [l  -  (Xi  +  -^)  -^] 


(A-93) 


'AP-5061-R 
Appendix  A 
Page  46 


In  Equation  (A>96) 


(A-97) 

(A-98) 

(A-99) 
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Roll 

Equation  (A-91)  may  be  written: 


[•*  +  ♦  -  +  V 


11 

K 


n 


(A- 100) 


In  Equation  (A* 100) : 


“0  “ 


a  (-gj-)  g 


D 

.D 

D 


(I  - 1 1  .to  ej)[  1  -  (x,  +  IP 

(A-101)D 

w 

D 
0 
n 
0 
D 
[ 


T  (0^ 

^  ,  JfxTT 


(A-102)' 


(A-IOS) 


AP-5O6I-R 
end  lx  A 
.^ge  48 


^iPesMfch  Manubciuring  Diviiion 

OmuIi,  Atimm 


APPENDIX  B 


FORMULATIGM  OF  GEM- SIX  DEGREE 
FREEDOM- BODY  MOTION  EQUATIONS 
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APPElfDIX  B 

FOKMULATION  OF  6Qi  SIX  DEGREE  FREEDOM  BODY  M0T1OH8 

I.O  GENERAL  EQUATIONS  OF  VEHICLE  MOTION  WITH  SIX  DEGREES  OF 
FREEDOM 

1,1  Velocity 

Consider  a  body  with  mass  (m),  body  fixed  axes  xyz. 

Inertial  system  XYZ  and  rotating  with  angular  velocity 
Thomson  (Ref.  2)  shows  that  the  absolute  velocity  of  any  point 
on  the  body  Is: 

^1  "  ^REL  +  CD  X  (B-l) 

where 

Vo  *  velocity  of  origin  of  xyz  relative  to  XYZ 
^REL  *  velocity  of  point  relative  to  xy* 

CD  X  r  *■  velocity  of  point  due  to  the  rotation  of  body 
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Since  xyz  is  a  body-fixed  axes.  It  Is  apparent  that  thcra 
can  be  no  velocity  of  the  body  relative  to  xyz.  Therefore, 
Equation  (B-1)  becamea: 


=  ^o  +  £0  X  (B-2) 

1.2  Momentum 

Since  linear  momentum  is  the  product  of  the  mass  of  the 
body  and  Its  Instantaneous  velocity  and  angular  momentum  is 
defined  as  the  moment  of  the  linear  momentum.  It  can  be  shown 
that  the  angular  momentum  of  the  body  Is  given  by: 

“  I  ^1  *  ®1 


=  S  r^  X  mj^  C^o  +  0)  X  r^^) 


=  Z  r^  X  mj^  (co  X  r^)  -  Vo  X  Z  m^^  (B-3) 


Integrating  over  the  body  and  noting  that  If  the  origin 
of  xyz  coincides  with  the  center  of  mass 


Z  r^  iDj  -  0  <B-^^) 
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Therefore,  the  angular  moiDentun)  written  as  components  in  xys 
become: 


Eo  =  h^  T  +  h  3  +  h^  ^ 


(B-5) 


Defining  the  moments  of  Inertia  about  the  xyz  axes 


as: 


ix  =  (y*  +  2* ) 

ly  =  ^ (x®  +  2*)  dm  ^ 

J  (x*  +  y*) 


(B-6) 


and  the  products  of  Inertia  as: 


1x2 


y* 


^  xy  dm 


X2 


Jyrd, 


(B-7) 


It  can  be  shown  that  the  angular  momentum  components  about 
XJ2  become: 
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h  S  I(JD  -  l  o)  -  lo) 

\  XX  xy  y  xz  z 

Sr  •  Vy  ■  •  V*".  } 


(B-8) 


Assuming  that  the  xyz  axes  are  the  principal  axes  (those 
axes  about  which  the  products  of  inertia  go  to  zero)  the 
angular  momentum  becomes:  * 


lo  .  1,01,1  +  lyOlyj  + 


(»-9) 


1.3  Moments 

Thomson  shows  that  the  moment  about  the  mess  center  of 
a  body  is  equal  to  the  time  derivative  of  the  moment  of 
momentum  about  this  point. 

Rq  =  Ho  -f  (0  X  Ho 


=  (l^T  +  hy3  +  h^H)  4  m  X  Ho  (•-10) 

The  cross  product  w  x  Ho  Is  shown  to  be  the  rotation 
of  vectors  h^^T,  hyj#  «nd  h^H  due  to  oj^,  cOy  and 
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or 


<0  X  Ho 

+  (<I>A 

Equation  (B-lO)  then  becooes: 

" 

Ro 

=  +  Hyl  + 

where 

•s. 

.  h,  +  «yh,  - 

“e'V 

• 

'  *V  ^  “e*^  ■ 

(D_h 

X  S 

► 

'  ^E  +  “eN  ■ 

“y'Sc 

4i 

Equation  (B-8)  substituted  into  Equation  (B-13) 
products  of  inertia  yields: 

H,  -  lyiy  +  VV”,  -  I 

-  Ve  *  va  -  Wv 


(B-11) 


(B-12) 


(»-13) 


neglecting  the 


(B-U) 
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1.4  Force! 

Thomson  also  denotes  the  force  acting  on  the  nest  at: 


?  -  FT  +  FT  +  ft 


(»-15) 


With  the  xyz  axes  system  being  body  fixed  and  rotating 
with  the  mass,  the  force  components  In  the  x,  y  and  z  directions 
may  be  determined  from: 


F 

II 

a 

+ 

u  X  mv 

These  components 

are 

!  seen  to 

be; 

F 

X 

• 

a 

+ 

v*<o„  - 

S  y 

V* 

F 

• 

=  ®(Vy 

+ 

- 

V  m 

y 

X  z 

Z  X 

F 

a  m(v^ 

+ 

V  m  - 

V  0) 

z 

y  * 

X  y 

(B-16) 


(B-17) 


1.5  Recapitulation  of  Motion  Equations 

The  general  equations  of  motion,  written  In  GEM  tamlnology 
become: 
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2.0  PARTICULAR  "SIX-DEGREE  FREEDOM"  EQUATIONS 

The  terms  comprising  the  six  equations  will  be  given  as 
"full  range"  parameters,  l.e.,  they  will  not  be  perturbated 
about  a  steady-state  point,  as  were  the  equations  for  pitch, 
roll  and  heave  In  Appendix  A.  The  terms  will  also  be  normalised 
In  such  a  manner  as  to  give  dimensionless  equations. 
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Let  the  vehicle  have  a  heading  V'  relative  to  XYZ 
Inertial  syateis  and  let  there  be  a  wind  at  an  angle  ^ 
with  the  Inertial  system.  The  relative  wind  can  then 
be  expressed  by: 


and 

Uy^  *  Uy  +  Sin  O  - 
Expanding  the  sine  and  cosine  terms  and  letting 
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FIGURE  B-1 
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=  b  (cos  3  cos  i/  +  sin  p  sin  V')| 

=  •  [^y  +  i  \  P  CO*  f  ■-  COB  p  sin  V')] 


(B-21) 

(B-22) 


The  relative  wind  will  be  used  primarily  as  a  drag 
force  and  will  therefore  appear  in  the  equation  as 
and  As  the  square  destroys  the  sign  of  the  terns, 

they  will  henceforth  be  written  as  follows: 


+  '7^(cospcosV<  +  sin^sinV')  (cos^cosV'  +  sin^sinV')j  (B-23) 

■*' a  ■  cospsinV')j|TJy  +  -  ^^^(sinpcosV'  -  cospslnv|^(B-24) 


The  parasite  drag  can  then  be  defined  as: 

=  "^D^B 

drag  coefficient  (c^)  here  is  related  to 
base  area  rather  than  frontal  area. 

1  -  fp  “*« 

and  from  Equation  (B-23)  yields: 
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[cos^coBif/  +  sin^8in'^')j^jj  +  V^(co8pcotV'  +  •ln^8ln^^  (B-25) 


Positive  Pitch  Angle 

=  -  Vb  ® 

It  Is  noted  that  In  the  steady -state,  Vb  very 
nearly  equal  to  the  vehicle  weight  (W^).  It  Is  therefore 
assumed  that  the  force  due  to  a  positive  pitch  angle  (nose 
up): 


-Wc  e 


since  6  Is  small,  sin  6=6  (radians) 
Let  0  “  ^  ® 


The  force  due  to  positive  pitch  angle 


becomes: 


Momentum  Drai 


(B-26) 


=  -Qf^xP 
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whar* 


=  total  volume  flow  to  fana 


and 


\  (cosgcosV'  +  sln^sinV'J 

fjj  =  \  (co8^co8V'  +  8in^8inV')j  (B-27) 

The  sucimatlon  of  forces  in  the  x  direction  equated 
to  the  X  component  of  Equation  (B-18)  yields  ; 


^x  +  ez  -  =  Wp  (^p  +  ^g) 


‘^dP^b** 


5-^ —  l^x  \(co8PcosV'  +  8ln^8inV')j  +  V^(co8Pco8^  +  sin^slnyr^ 

2hWc  -  r  1 

— ^  e  -  Q^pb  1  {cos^cosf  +  8in^8in|(')J 


Let  X 


where 


W  ' 


^x 

T-> 


u 

_2 


and  neglect  ez  because  of  relative  magnitude. 
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dT) 


IT  = 


•dps^b^ 

“■SW_ 


■'■  ^^^(cosficos’yi-  +  aln^sinV/)|^j^  +  ^^{cospcosV'  +  slnpaln^'^ 


2hg  ^ 


QfP8 


[|^x  ■'■  (cos^cos^  +  sln^siny/j 


(B-28) 


yr 


=  a‘ 


2.2  Side  Velocity 
Parasite  Drag 

The  relative  wind  Is  given  In  Equation  (B-24)  as 

l^y  +  j  \  (sln^cosV'  -  co8^slnV')j ^TJy  ■*■  a  \  (slngcos^'  -  cospsln^)| 

The  parasite  drag  then  becomes: 

£  .  U 

y  2a  yr 


|r  .  I  (.inpco., 


-  c 


os^slnV')  IK  +  j  (slnpcosV'  - 


:o8^8lnV')j 

(B-29) 
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It  Is  noted  that  this  term  contains  the  saioe  drag 
coefficient  as  the  equivalent  term  In  the  foirward  velocity 
equation.  It  Is  therefore  modified  by  the  ratio  of  planforn 
length  to  width  (b/a). 

Positive  Roll 


=  Vb  =  '^c  ^ 


(B-30) 


where 


0  =  5E  ^ 


Momentum  Drat 


fy  =  -  QfpUy  =  •  QfP®^^y  +  i  \  (slngcosV'  -  cos^slnV'^  (B-31) 

The  summation  of  forces  In  the  y  direction  equated  to 
the  y  component  of  Equation  (B-18)  yields: 


+  V*  -  =  -  QfP«^y  i  (sln^cosV'  -  cospslniff^ 


_£  J  — iTIy  +  j  V^^(sln&co8^  -  cosfsslnV')j  [iJy  +  j  ^^^(slnpcos^  -  cospaln^] 
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Let  y  = 

at*' 

where 

U 

II 

a*. 

\  -  f 

and  neglect 

<bz  because 

of  relative  magnitude. 

dTi 

- 


QfP8 


^TJy  +  ^  (sln^cos^  -  cospsinV')J  + 


-  [Uy  +  i  V^(sln^cosV'  -  co8psin^)J^y  +  j  \{s±T)^coBi/ 


co8p8in^J-j 

(B-32) 


2.3  X«H 

The  yaw  component  of  Equation  (B-I9)  Is: 

”z  =  h  *  *  ®  * 

The  external  moments  that  constitute  the  left  side  of  the 
above  equation  will  Include  the  differential  thrust  on  the 
propellers  and  aerodynamic  forces  due  to  side  velocity  (relative 
wind). 
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Differential  Thrust 


Let  be  the  moioent  arm  of  the  propellers  and  the 
moment  about  z  due  to  propeller  thrust  becomes: 

-  "c*p  <’p  -  ’,)  ‘*-33) 


It  Is  reasonable  to  assume  that  any  properly  designed 
vehicle  will  have  a  certain  amount  of  yaw  stability.  This 
will  arise  from  aerodynamic  forces  due  to  a  side  velocity 
(side-slip)  which  produces  a  moment  opposite  to  the  angle 
of  yaw.  It  is  quite  obvious  that  this  "weather  vane" 
effect  must  be  built  In  the  vehicle,  lest  the  operator 
be  continually  "riding"  the  yaw  control. 

If  this  force  Is  defined  as  a  drag  force  due  to.  a 
side  velocity  It  Is  seen  that: 

"z  =  'D  <5  PX'yr'^B  ?  (»-3'*) 

where 

b^  =  "moment  arm"  of  relative  wind  (aft 
of  vehicle  c.g.) 

and 

U  =  relative  wind  given  In 

JT 

Equation  (B-22) 


AP-5061-R 
Appendix  B 
Page  17 


^iReseardi  Manu^uring  Division 

HumlM,  Jiimm 


These  external  noments  equated  to  Equation  (B-19) 
yield: 


Htt 


a_gW 
p  c 


(’p 


+  a  ^|,(sinPcosV/  -  cospsin^)j  iJy  +  j  ^^(8in^cosV'  -  cospsin’^^  (B-35j 

2.4  Pitch 

Equation  (A-96)  of  Appendix  A  gives  the  normalized, 
linearized  equation  for  pitching  motion  due  to  variations 
In  fan  flow  and  wave  disturbance.  Neglecting,  for  the  present, 
the  wave  disturbances,  the  equation  is: 


p'  +  2Cg  .  +  co/j  e  = 

From  which 

It?  =  ■  <^e  H  -  0 

■  ““0  ^  "  ^fa  ■  ”^fs 

O0fi. 
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The  external  moniehts  indicated  by  My  in  Equation  (B-I9) 
and  which  will  apply  also  to  the  above  equation  will  be  shown 
to  include  moments  arising  from  aerodynamic  forces,  propulsive 
forces  and  surface  displacements. 

Up  to  this  point  in  the  formulation  it  has  been  assumed 
that  the  surface  over  which  the  vehicle  was  "flying"  was  rigid, 
that  is,  there  has  been  no  allowance  for  the  surface  to  deflect 
due  to  a  change  in  surface  pressure  as  will  occur  whei.  a  vehicle 
"flies"  over  water.  This  characteristic,  which  will  be 
analytically  described  in  detail  in  Appendix  C,  is  a  necessary 
parameter  in  a  complete  formulation. 

It  is  shown  in  Appendix  C,  Equation  (C-52)  that  the  mean 
slope  of  the  water  surface  as  a  function  of  forward  velocity  Is: 

-  cos  a) 

where 

It  can  be  seen  that  as  a  vehicle  "flies"  over  the  water 
surface  at  some  forward  velocity  (U^)  it  will  tend  to  follow 
the  shape  of  the  surface  deflection,  that  is,  as  forward  speed 


AP-506I-R 
Appendix  B 
Page  19 


^iGescarch  Manulicturing  Division 


Mmwi,  4Htm* 


Increases,  the  vehicle  will  tend  to  nose  up  due  to  the  formation 
of  an  uphill  slope  beneath  It.  It  will  continue  to  pitch  up 
until  the  "hump"  speed*  Is  reached  at  which  time  the  vehicle 
will  return  to  a  near  horizontal  altitude. 


As  this  water  slope  (e^)  can  be  thought  of  as  a  disturbance. 
It  can  be  Included  In  Equation  (B-35).  Normalized  to  Equation 
(B-36)  the  disturbance  becomes: 

“‘e  ”c  ri2  «lna  ,»lna  ...  0  /. 

Spjlgh  -5—  (“a  -  “•  “2J 

In  addition  to  the  above  description  of  water  surface 
behavior.  It  will  be  well  at  this  point  to  Include  the  term 
which  describes  the  transient  behavior  of  the  surface  due  to 
pitching  motions.  This  term,  as  fully  described  in  Appendix  C, 

Is  seen  to  modify  the  term  of  Equation  (B-36)  by 

-  ^)  <“'9 


*"Hump  speed"  Is  defined  as  the  speed  at  which  the  slope 
of  the  water  Is  a  maxlmuD. 
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AHimt 


wh«r« 


16  I 


«  p7^ 


2  sin  a 


cos  a) 


1+0)® 


16  I. 


(B-38) 


Pitching  Moment  due  to  Propeller  Thrust 


”e  -  -S 


Zp  =  moment  arm  of  propellers 


Normalized  to  Equation  (B-36)  the  above  term 
becomes: 


-Iei-  (’p*’.) 


(B-39) 


Pitching  Moment  due  to  Aerodynamic  Lift 

If  the  aerodynamic  lift  does  not  act  through  the 
vehicle  c.g.  and  It  can  be  shown  that  It  does  not,  ttere 
will  occur  a  pitching  moment  due  to  forward  velocity. 


Defining  the  aerodynamic  lift,  as 

h  ■  'l  <5  P  "*)  *» 
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For  the  pitch  node,  only  the  forward  velocity  will  I 

be  considered.  — 

With  the  relative  wind  from  Equation  (B-23)  and  1 

normalizing  to  Equation  (B-36)>  the  aerodynamic  pitching  | 

moment  becomes: 

V^^(cospcosV'  +  sln68ln:^)|  +  ^^(cospcos^  +  slnpslnV')|  (B-40)^ 

Appendix  C  shows  the  variation  of  base  pressure  due  to  a 

forward  and  side  velocity.  This  term  which  produces  an 
Increase  In  base  pressure  In  the  forward  compartments  and  | 

a  decrease  In  pressure  in  the  aft  compartments,  produces 
a  nose  up  (+d)  moment.  ( 


The  variation  in  base  pressure  depends  upon  the 
relative  wind  as  do  the  aerodynamic  forces,.  Normaliz¬ 
ing  to  Equation  (B-36)  this  moment  becomes: 


pgAgb*  (.29«) 


4hIy4(a+b+2aT 


~|0x  +  V^icoa^coBtj/  +  sinbslnV')|^jj  +  V^(co8bco8V'  +  8in^slnl^^  || 

'"‘I 
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The  summation  of  the  above  outlined  pitching  moments  whan 
substituted  Into  Equation  (8-36)  yields: 


It®  =  '•  ^  -  (1  -  Ig)  £i)*g  ? 


-  OJ 


r 


35, 


<’£,  -  •  ’f,  +  V) 


2  sin  a  /Sin  a  cnn  '  ^o****c 

“  e  p7?  -  cos  - '  -e—  - 


+  4(a+^lar^'^  I  ^  ^  os&c  08  v>+8  ln^slnV'  )|  ( co8pcos»+s  tnps  tn^j 


[I 

fl 


(8-42 


2.5  Roll 


Equation  (A-lOO)  of  Appendix  A  shows  the  normalized,  linearized 
equation  for  rolling  motions  as  a  function  of  fan  flow  variations. 
Rewriting  this  equation: 


d*f 

Ht5 


§  - 


CD 


*  ? 


0 


+  "%  ^  <’£i  *  'fs  -  "£.  -  '£.> 


(»-«) 
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The  external  rolling  moments  indicated  in  Equation  (B-I9} 
are  assumed  to  arise  only  from  aerodynamic  forces. 

The  change  in  water  surface  slope  is  neglected  because  it 
is  not  anticipated  the  vehicle  will  travel  at  high  speeds  in 
the  side  direction.  It  is  further  assumed  that  "side  slipping" 
at  high  forward  speeds  will  occur  above  the  hump  speed.  The 
transient  behavior  of  the  surface  due  to  rolling  moticms  will 
be  included  however. 

As  there  might  exist  crosswinds  of  appreciable  magnitude 
the  aerodynamic  moment  due  to  relative  wind  will  be  included 
as  well  as  the  moments  due  to  an  increase  in  base  pressure  et 
forward  velocity. 


As  in  pitch  the  term  which  describes  the  transient  behavior 
can  be  shown  to  modify  the  term  of  Equation  (B-43)  by: 


(1  - 


where 


16  I 


0) 


X 


<t> 


(1  +  €  -  ^  COS  a) 


■«) 


1  +  a)‘ 


s 


16  I 


(B-H) 


h 
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n 

fi 

n 

n 

0 


vioa 

Wmit.  j<i<— 

This  tern  Is  more  fully  described  In  Appendix  C. 

With  the  relative  wind  from  Equation  (B-24)  and  normalizing 
to  Equation  (8-43)^  the  aerodynamic  rolling  moment  becomes: 


- 


c.pgA^®by  I  V  If  b 

- lETT? - a  +  *  V^^Csin^coaV'-cospelni^j 


(B-45) 


As  shown  In  Appendix  C  the  variation  of  base  pressure  due 
to  velocity  produces  a  rolling  moment  as: 


-  pgA.a®b(0.365a)  I  IL  k  1 

'  4ht~4(«4.bt8b?i j~ Py  *  •  +  i  'l',(.tbeco.»-COiPstn»)J 


Equation  (B-43)  then  becomes: 


aP  -  -  ^  3^  -  «£,  -  ’f.) 


pgAg.»b 

Tir^ 


0.365a 

4(a+b+2bT 


— j  |TIy+^^(slnpcosV'-cos^8in^)|^y+j  l^^(sln^cosilrTCOspsln^) 


(•-♦7) 
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It  is  noted  that  the  products  fe  and  were  neglected 
because  of  their  relative  magnitude  In  the  pitch  and  roll 
equations. 

2.6  Heave 

The  normalized 4  linearized  equation  of  motion  In  the  heave 
mode  from  Equation  (A-92)  (neglecting  wave  disturbance)  la: 

^  -  2^z  ^  ‘  “*z^  -  “*z  ^^fx  +  ^fa  ^fa  ^£4^ 


As  shown  In  Appendix  C,  the  frequency  term  In  the  above 
equation  will  be  modified  by  a  term  that  shows  the  transient 
behavior  of  the  water  svurface  due  to  heaving  motions. 


The  term  becomes  (1  -  T^)  a)*^l 


where 


^  (1  +  €  -  ^  coa  a) 

- jg - - 


(B-W) 
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As  stated  previously ^  the  lift  due  to  the  trapped  pressure 
Is  augmented  by  an  aerodynamic  lift  arising  with  velocity.  It 
Is  assumed  that  the  lift  due  to  side  velocity  Is  small  and  may 
be  neglected.  Hence,  only  the  forward  velocity  will  be  considered. 
It  Is  further  assumed  that  the  variation  of  base  pressure  due  to 
forward  velocity  Is  nearly  balanced  fore  and  aft  and  may  be 
neglected, 


The  aerodynamic  lift  term  Is: 


"5W  h"  ^xr 


or  with  from  Equation  (B-23) 


l^x  \  (cospcosV'  +  slnpslnV')j^j^  +  \  (cospcosf  +  sln^sln^)j 


(B-II9) 


The  equation  for  heave  motion  then  becomes: 
(with  additional  terms  from  Equation  (B-18)  ) 
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dZ 


^^{eoB^coaf  +  8ln^8ln^)j  +  V^{coa^cqaf  +  8in3sln^^ 
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APPENDIX  C 


SURFACE  DEFLECTIONS  AND  BASE  PRESSURE 
VARUTIONS  DUE  TO  FORWARD  VELOCITY 
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MUmm,  Aritmt 

SURFACE  DEFLECTIONS  AND  BASE 
PRESSURE  VARIATIONS  DUE  TO  FORWARD  VELOCITY 


1.0  Surface  Deflections 

To  more  fully  describe  the  behavior  of  a  vehicle  hovering 
and  maneuvering  over  water  it  was  necessary  to  investigate  the 
behavior  of  the  water  surface  when  subjected  to  a  uniform 
pressure  over  the  surface.  The  following  then,  presents  a  two 
dimensional  problem  solution  of  standing  waves  in  a  stream  of 
infinite  depth  due  to  a  uniform  pressure  distribution  over  a 
finite  length  of  the  surface,. 

The  general  solution  of  the  problem  is  given  by  J.  J.  Stoker* 
in  his  book  Water  Waves  in  the  form  of  an  integral  of  an 
analytic  function. 


List  of  Symbols  (used  in  this  appendix  only) 


X 

y 

t 

n(x,t) 

u 

Po 

*Reference  3 


horizontal  space  coordinate 
vertical  space  coordinate 
time 

surface  displacement  (from  y  =  0) 
free  stream  velocity 
pressure  on  the  surface 
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List  of  Symbols  (Cont'd) 

P  -  water  density 
g  -  gravity 
0  -  velocity  potential 


FIGURE  C-1,  POSITIVE  SENSE  OF  COORDINATES 

In  brief,  the  problem  is  formulated  as  follows:  (Sea 
Chapter  7  of  Reference  3)- 

The  velocity  potential  roust  satisfy  the  LaPlace  Equation, 


-  <t>  +  (#1  ■  =  0 

V  ^  ^xx  ^yy 


(C-1) 
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Mmiws,  Arimum 

The  value  of  obtained  from  Equations  (C-3)  and  (C-4) 

~ satisfies  the  boundary  conditions  mentioned  as  well  as  the 
condition  that  the  surface  pressure  is  po  over  the  region  froo 
-a  to  -fa. 


FIGURE  C-2 


Stoker  states  that  Equation  (C-3)  must  be  Integrated  from 
+  i  <*o  (positive  imaginary  axis)  to  any  point  (Z)  in  the  lower 
half  plane  (negative  y).*  Interest  here  is  the  surface  shape 
which  will  be  obtained  by  evaluating  rj  from  the  first  of 
Equation  (C-2).  Noting  that  and  r|  are  independent  of  time 
we  obtain 


(C-5) 


*The  path  of  integration  must  inclrcle  the  origin  and  the  points 
(-a,  o)  and  (-m,  o)  in  a  counterclockwise  fashion. 
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MmMX.  ilfiMM 


Since  the  equation  for  q  contains  not  <i>  but  Its  derivative 
with  respect  to  x.  It  will  be  convenient  to  differentiate  £(Z) 
with  respect  to  x. 

Equation  (C-4)  then  yields 

‘  ♦x 

It  Is  noted  that  partial  differentiation  of  a  function 
of  Z  with  respect  to  x  Is  the  sane  as  total  differentiation 
with  respect  to  Z  because: 

Differentiation  of  Equation  (C-3)  yields: 
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is  z 

^og  ^  0 

TJ3^  ®  J 

loO 


logffi  it 


(C-9) 


Integrating  the  first  term  by  parts  gives: 


FTi) 


dt 


(C-IO) 

Substituting  Equation  (C-lO)  Into  (C-9)  yields: 


too 


(C-ll) 
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fkmUMt  Afimm 

It  Is  seen  that  the  latter  two  terms  will  cancel.  The  tern 
remaining  to  be  Integrated  does  not  yield  elementary  functions 
but  rather  yields  an  Infinite  power  series  upon  successive 
Integration  by  parts.  Therefore,  to  obtain  numerical  solutions. 
Equation  (C-11)  will  be  integrated  numerically.  Letting  b  denote 
a  particular  point  on  x.  Figure  C-3  shows  the  four  regions  of 
Interest. 

Also,  define: 


a 


T)  = 


1 

a 


€  =  T 


(C-12) 

(C-13) 

(C-U) 

(C-15) 
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Define  five  paths  of  Integration: 

(1)  Let  t  =  ly,  dt  =  Idy 

>  y  ^  a 

(2)  Let  t  =  x  +  la,  dt  =  dx 

0  ^  X  ^  b 

(3)  Let  t  =  b  +  ly,  dt  =  idy 

a  ^  y  ^  0 

(4)  Let  t  =  X  +  la«  dt  »  dx 

0  ^  X  ^  -b 

(5)  Let  t  =  -b  +  iy,  dt  =  idy 

^  1  y  1  0 

Figure  C-3  shows  the  possible  paths  of  integration. 

Paths  (1),  (4)  and  (5)  will  determine  (f>^  for  -b,  p  >  1. 

Paths  (1),  (2),  (3),  (4)  and  (5)  will  determine  0^^  for  -b* 

P  <  1;  +b,  t  y  1  and  +b,  p  <  1  provided  the  dashed  (closed) 
paths  shown  are  added  to  (1),  (2),  (3),  (4),  and  (5). 

.  r 
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In  view  of  Equations  (C-3)  snd  (C-Il)  the  function  to  be 
Integrated  la: 


r  £5 

X  ~  po 


1  z 


z 

a 

I 


I  - 1 


d  ( 


(C-16) 


The  five  Integrals  now  become: 


Path  (1) 


'(1) 


'  ■  u. 


.-an 


mjVi  -  Tn-TT^ 


idq 


i  r  O’. 


=  + 


I  O’. 
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(C-17) 


Path  (2) 


-  -if  (t-^4^  ■  r^fn) 

0  ^  ' 

=  ■  4  ®"°  5  8ime)4f5^^  -  (l+ijrn:^ 


K 


co8ae+  i  sime) 


€  +  1 


/  1  _  1 
he  -  1)2  Ti  -  rrTT 
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12) 


-a 


S  ^€-1)^  -  '(rniW  ■  (€+1)«T^' 


C08a€ 


-i 


e  dc 


^  (tApti  ■  tAt^  ■  (€+i)4+i)““‘j  '“'j 


(C.18) 


Path  (3) 


13) 


=  ^  gla(B+in) 


lob 


S  “-“"(tH 


-1-iTi 

TTW 


B+l-ln 


^  dn 
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Anmm 


'(4) 


■m 


e  -  1 
i(e+1)*+1 


ITiTF 


stmc  -  ^ 


T€-r)»TI  - 


+  i 


0  L 


€  -  1 

lTe-l)2-,T 


(C-20) 


de , 


Path  (5) 


'(5) 


w 


,la(-P+ln) 


/  1  1  \ 
I -e-f-iTj-i  -s+lTi+I j 


IdTJ 


-iaP 


e 


-ari 


(e+i)i+'n2  • 


dTi 
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1 

It 


cosaP 


A 


6-1 


6  +  1  \ 

4(6+1F1  • 


s'(. 

o  \ 


+  8ina6  \  (  -  n5'Zf?4.nr)  ®  ^ 


i+l6+l)*) 


] 


-1)^  ■ 


6  +  11 
4(6+1)“  J  *  ‘’t 


cosae 


K’ 


Ti2+(e-ij2  n® 


e 


(C-21) 


In  order  to  define  the  closed  paths  (dashed  paths  of  Flg- 
ure  C-3)  needed  to  supplement  the  five  preceding  Integrals 
when  integrating  to  one  cf  the  four  regions  of  interest  we  will 
make  use  of  Cauchy's  Integral  Formula. 

Caucl^'s  Integral  Formula  states: 

If  f(Z)  is  regular  in  a  region,  then  the  formula 
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'(z)  .  jIi  ^  ^  dc 

Is  valid  for  every  simple,  closed,  positively  oriented 
path  for  every  point  Z  in  its  interior. 

or 

^  dC  =  2Tif(Z)  (C-22) 

For  example: 

let  f(C)  = 

C  =  t,  dC  =  dt 

Z  *  -1 

therefore 

■  -i  §  -  -Sty'll  I-*-**) 

=  -2e"^®  =  -2  (cosa  -  i  slna ) 

where  -  ^  is  an  arbitrary  multiplier. 
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In  Equations  (C-l?)  through  (C-21)  let: 

-  I  ^  ({eUj'il  ■  ) 

•  i  ^  (tAt4i  '  T£+1)»1iJ 


'♦  -  i  J  ■  Ie+l)»+l) 


coaa€  de 


(C-23) 


n^+(e 


V  <"> 


New  with  the  five,  paths  of  Integration,  the  Integrals 
themselves  and  the  closed  paths  defined  we  may  define  the 
function  for  the  four  regions  of  Interest. 
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For  the  region  Z  =  -b,  S  1  (Figure  C-3B) 


4  =  e-a^ 

X 


Ii  +  (Fi  -  F4)  e-a  +  i  (Fa  +  Fs)  e-» 


-  Gi  co■a^  -  G3  slna6  +  i  (Gi  sitiap  -  Gs  cotaS) 


if  =  coaa?  +  1  stnafi 


=  Ix  coaaP  +  e“®'  (Fi  -  F4)  coaafi  -  e-a  (p^  +  Fa)  •imp 


Ga  comp  simp  "  Gx  cos^aP 


+  Ga  comp  simp  -  Gx.8in*aP 


?  =  Ix  comp  -  Gx  +  (Fx  -  F4)  e"*comp 


(Fa  +.  Fa)  simP 


(C-24) 


For  the  region  Z  =  -b,  P  <  1  (Figure  C-3A) 


Equation  (C-24)  will  apply  for  this  region  provided  the 
following  tern  fron  Cauchy's  integral  is  added: 
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Am 


<t>^  =  Ii  cosQ^  -  Gi  +  (Fi  -  F4)e”*  cosap  -  (Fa  +  F3)e 


The  additive  term  due  to  Cauchy's  integral 


^-la  ^  .  2,-ia(H-?)  _  .  ^  (l^g) 


i  2  Sim  (U^) 


Therefore  > 

=  Ii  coao^  -  Gi  +  (Fi  -  F4)e”°'’  cosa^  -  (Fa  +  Fs)e“*  sina^ 

.*  2  cosQ  cosa^  +  2  aim  sim^  (C*26} 

For  the  region  Z  =  +b,  ^  >  1  (Figure  C-3D) 

The  same  open  path  as  the  previous  region  will  apply  here 
if  the  following  closed  path  is  added: 

(-  2.-^^  +  (2.-^  ,‘»)  =  *  »  .i«  ,lnaP 

therefore 2 

/ 

«  Ii  cosQ^  -  Gi  +  (Fi  -  F4)e"'®  cosa^  -  (Fa  +  F3)e"*  sinaS 


i  f  4  sina  sinaP 
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UtmbkAitmm 

As  stated  previously^  the  integral  given  In  Equation  C'll 
and  hence  the  integrals  occurring  in  Equations  (C-24)  thru  (C-27) 
will  not  yield  an  elementary  function  but  rather  an  infinite  power 
series  upon  successive  integration  by  parts. 

In  order  to  numerically  integrate  the  functions  of  interests 
Simpson's  Rule  was  applied  to  the  integrals  and  a  digital  computer 
program  developed  to  evaluate  the  functions. 

It  is  obvious  that  the  expressions  derived  in  the  foregoing 
section  are  far  too  complex  to  be  of  any  use  other  than  as  a 
comparison  for  the  approximate  formulation  in  the  following 
section.  Therefore^  for  purposes  of  comparison^  several  wave 
shapes  were  determined  utilizing  the  above  foinmulation  and  the 
Digital  Computer.  These  surface  shapes  are  shown  graphically 
in  Figure  (C-4). 
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1.1  simplification  of  Wave  Shape 

To  obtain  a  more  convenient  expression  for  the  wave  shapes 
the  following  equation  was  developed  which  Is  seen  to  very  closely 
approximate  the  shapes  shown  In  Figure  C-4. 

(C-28) 

where 


€  here  Is  a  function  of  a  (see  Equation  C-50) 


Zw 


=  -f 

PS  L 


1  +  €  -  2  cosa  (1 


-f)] 


FIGURE  C-5 
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Equation  (C-28)  represents  a  traveling  wave  of  fixed 
shape  moving  with  the  craft.  Note  the  change  In  coordlnatea 
In  Figure  C-5* 

The  average  vertical  displacement  of  the  surface  is 
obtained  by  Integrating  Equation  (C-28)  from  x  =  -b  to  x  -  b 
and  dividing  by  2b. 


^  1+€-2  cosa  (1  -  |)  dx 

(1  +  e)  2b  +  2  I  sina  (1  -  f)  . 


Zo  =  ^  (1  +  €  -  ^  (C-29) 

The  general  slope  of  the  water  surface  will  be  obtained 
by  finding  the  slope  of  the  best  straight  line  approximation 
to  Equation  (C-28) .  This  Is  done  by  making  the  Integral  of 
the  square  of  the  difference  a  minimum. 


P8 
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?o 

pg 


(C-30) 


-  Z^)*  dx  (C-31) 

B  will  be  obtained  if: 

Z^)  dx  -  0  (C-32) 

Z^)  X  dx  =  0  (C-33) 

•)  and  (C-30)^  Equations  (C-32)  and 

-  €)  +  Bx  +  2  cosa  (1  -  ^  )"|  dx  «  0  (C-3^) 
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dl 

3F 


s::  [' 


(a  -  1  -  6)x  +  bx*  +  2x  C08Q  (l  -  g)  ]  dx  -  0  (C-35, 


.[ 


dx  =  b  -  (-b)  =  2b 


X  dx  =  i  (b®  -  b*)  =  0 


d*  .  1  b>-  (-b) 


2  vs 


(C-36) 


(C-37) 


(C-38) 


cosa  (1  -  5)  dx  =  +  ^  sin  2a 


a 


(C-39) 


S-b  *  (1  -  f)  dx  =  I?  U  -  cos  2a  -  2  sin  2a)  (C-AO) 

Using  Equations  (C-36)  through  (C-40),  Equations  (C-^) 

•nd  (C-35)  become: 


2b  (A  -  1  -  €)  -  2  I  sin  2a 


or 


A  =  1  +  € 


sin  2a 


(C-Al) 
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fl 

0 


■  .U 


and  2a  -  a  tin  2a)  «  0 

or  . 

B  =  -  (1  -  coa  2q  -  Q  8in.2Q) 

B  =  _  ^  •isa.  (Sisa.  -  coaa)  (C-42) 

D  a  a 

To  obtain  the  maxlcaim  slope  Equation  (C-42)  Is  differentiated 
with  respect  to  a  and  equated  to  aero. 

This  yields: 

*  3  +  2  (1  -  a*)  (C-43) 

The  lowest  root  of  the  above  equation  Is: 
a  =  1.753.  =  0.561,  cos  a  -  -0.1812 

The  maximum  value  of  B  becomes 

®max  =  ■  F  ^0-561)  (0.561  +  0.1812) 

-  -  ^2®  (C-44) 
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Therefore,  the  maximum  slope  !•: 


dZ 


max 


Pa  B  - 

pg  max  pgb 


(C-45) 


Numerical  value  of  e  — . 


The  parameter  e  is  a  function  of  a  and  is  obtained  by 
making  the  approximate  expression  for  water  surface  given  In 
Equation  (C-28)  conform  as  closely  as  possible  to  the  exact 
value  given  In  Section  1. 


Let  the  exact  solution  be  defined  as: 


r)  -  (C.46) 

and  the  approximate  value  be  defined  as  In  Equation  (C-28) 

-  1  -  2  COS  (1  -  |)  +  €  (C.47) 

Define  now  the  sum  of  the  squares  of  the  variation  as; 

•  ®  (C-48) 


The  I's  represent  N  distinct,  uniformly  spaced  points  at 
which  and  t]  are  evaluated.  It  Is  seen  that  making  £  a  mini¬ 
mum  provides  the  best  possible  value  for  e:  This  is  done  by 
making  the  first  derlvating  sero. 
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Expanding  the  square  of  Equation  (C-48)  yields: 


L 


and 


ds  _ 

ar-  2 


-  Ij)  +  2«  .  Tij)  +  c» 


(0-49) 


or 


-  T)) 


(C-50) 


The  values  of  e  for  various  values  of  a  are  deterinlned 
using  results  of  the  exact  solution  from  the  Bendlx  Digital 
Computer.  These  values  are  shown  graphically  in  Figure  C-6, 
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a 


FIGURE  C-6 
€  VS.  a 


fl 

D 

D 

D 

B 

D 

D 

C 

[ 

[ 

[ 


t 

1 


The  mean  displacement  (A)  and  the  mean  slope  of  the  water 
surface  (B)  in  dimensionless  equations  now  becoae: 


A  ■  1  +  e  -  2  cosa 

bB  -  6  -  coBa) 


(c-5i: 

(C-58: 
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2.0  Effects  of  Relative  Wind  Upon  Base  Pressure 

To  further  describe  the  effects  of  forward  and  side 
velocity  (relative  wind)  upon  vehicle  characteristics,  the 
induced  pressure  around  the  vehicle  due  to  relative  wind 
will  be  Investigated. 

2.1  If,  in  Equation  (1)  (AP-5047-R),  an  aabient  pressure 
Is  Included  on  the  right  side  of  the  equation,  the  end 
result  can  be  seen  to  be  the  redefining  of  the  parameter  5, 
such  that 

,  .  (l*!  - 

P  Vj* 

From  Equation  (A-33),  the  following  steady-state  equation 
can  be  written  for  the  peripheral  jets,  in  general, 

7T,  \  -  2Kj^  (1  -  f  K  (Po  -  P.i)  (C-53) 

C 


Also  for  compartment  Jets  from  Equation  (A- 37) 


_ik  n 

q  \  Qlk 


■  Slk  ^^1  ’  ^Ik^ 


Ilk  ■  "'plk 


In  the  above  equations,  note  that 


^1  ^ 


-  i  « 

Ilk  2 
c  .  •  ®lk 

Sik  cr 


(C-54) 
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In  the  steady  state 


-  it* 

TPa  5i2 

TTg  Qgi 

-it* 

PTg  ^2 

5.  - 

TTS  ^41 

+  TTa  ^14 

+  TTa  Qes 

+  TTa  ^4 

+  'T.  54. 


(C-55) 


If  Equations  (C-55)  are  expanded  using  Equations  (C-53) 
(C-54)j  they  will  yield  a  steady-state  equation  for  each  base 
sector.  Assuming  continuity  (the  net  flow  out  of  a  sector  is 
zero)  when  the  ere  zero  and  also  a  increase  on  each 

pressure  when  not  zero  yields  the  following  set  of 

equations : 


Let  fi 


ai 


h  ^al 


(i  +  rwA. -f.  §^s. 

(1  +  f  j  ^  A,  -  f,  it*  S, 

(1  +  pj  A,  -  f,  it*  ^  if-  ■  ^11. 

(1+P4)S4-?4§J*A,  -f,§^s..p„ 


(C-56) 
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Letting  A 


ClA _ 

C^(l  +^1) 


_ — 

Ci(i  +  Ti) 


The  solution  of  written  in  matrix  form  is: 


B(l+A*-B*)  (1-A2-B*)  A(l+3*-A*)  •  2AB' 


(l-A*-32)  B(1+A2.B*) 


A(l+32-A2)' 


(1+3* -A*) 


A(  1+3* -A*  )  ( 1-A*  -3* )  3 ( 1+A* -B* ) 

2 AE  B  ( 1+A*  -3*  )  (1-A*  -?*  ). 


P-  /(l+FiT 

P,^/(1+Pi) 


where  ^dJ  »  determinant  of  coefficients. 


(C-57) 


c  c 

Assuming  that  acts  at  the  center  of  the  sector  2 

the  equations  for  pitching  and  rolling  moments  are  seen  to  be 


C 

c  j 


(-  Aj.  +  Aa  +  A*  -  A4) 


- — 

(1  +  2F.  §^)  '■ 


+  F.  +  P.  - 

as  as  a+J 


(C-58) 


AP-506I-R 
Appendix  C 
Page  33 


and 


diCMMfdi  Manu^uring  Divsion 

fkmtkiat  MHmm 


c 

— fl  (Aj.  +  A*  -  A#  -  A4) 


T 


iXS. 


(1  +  sf.  g^) 

2 .2  Assumed  Ambient  Variations  Due  to  Forward  Velocity 

The  assumed  induced  pressure  coefficients  are  as  followa: 


Due  to  forward  velocity  (+  Uj^) 


(C 


D 
E 
E 

-59)[] 


front 

+  0.35 

sides 

-  0.23  q* 

rear 

-  0.23  \ 

Due  to  side  velocity  (+  Uy) 

front 

-  0.23  qy 

rear 

-  0.23  qy 

port  side 

-  0.23  qy 

starboard  side 

+  0.5  qy 

(C-60)' 


0 

D 


(C-61) 


D 


where  j  P 

qy  -  2  P  Vy* 
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!*«■)■.  ditmm 


-  0*83  q_ 
111  111 


0.23 


f 

a/2 

_ 

+0.35 


-  0.23  q- 

111  111 


111  1  III  111  i  III 

y  "  0.23  qjj  7  +  0.5  q^ 


POSITIVE  FORWARD 
VELOCITY 


POSITIVE  SIDE 
VELOCITY 


FIGURE  C-7 

INDUCED  PRESSURE  COEFFICIENTS 


In  pitch: 


0.23  q. 


0.35  f  -  0.23  I 
«  +  b 


o.35«-  --A23b 
TTb 


<lx 


Equation  (C-58)  then  becomes 


am. 


T 


iXA- 


(1  +  27x 


+  2(0.23)  q*  +  2 


recalling  that  ^ 

C 
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r 


lltll 
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aiwuii,  AHimt 


AM. 


T 


Qt23-g>. 


(1  +  ^i) (*  +  b) 


(C-fe) 


In  roll: 


p  K  p  =  ~  0«2'^*  q 

as  a  +  b  y 


Equation  (C-59)  becomaa: 


AM. 


B” 


Cn 


(1  +  §t 


) 


2(0.23)  il„  - 


.St  -.  p.a^l  •  4 

a  +  b 


^  0r^65  ab 

-  T  „  .  2b  ^ 


]r7||7Fr)(r+b)  y 


(C-63) 
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Effect  of  Water  Dvnamtca  Upon  Body  Motj 


As  a  vehicle  pitches^  rolls,  or  heaves  there  are  restor¬ 
ing  moments  and  forces  generated.  It  Is  assumed  that  In  the 
pitching  mode  this  restoring  moment  Is  given  by  the  following: 
(See  Figure  C-8) 


FIGURE  C>8 


(I)  (I) 

-  -  ^  AP  Ag  b  (C-64) 

Denote  this  restoring  moment  as  then 

-  -  -  I  AP  Ag  b  (C-65) 

9  above  Is  measured  from  a  horizontal  surface.  If  the  vehicle 
Is  over  water  this  angle  will  not  be  simply  9  but: 
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A0  -  0  -  e_ 


Where  ■ 
w 


4AP 


Equation  (C-65)  becomea: 


-  "9  (  «  -  B^)*  -f 


or  Ag  b  AP 


Kg0 


1  + 


Tire 


(C-66) 


In  elementary  terms  the  equation  of  motion  In  pitch  may  be 
written: 


I  My  -  -  I  ap  Ag  b  +  My'  -  ly  9 
where  are  external  momenta 


(c-er) 


and  substituting  Equation  (C-66)  Into  Equation  (C-67)  yields: 


•  # 

0  + 


— ~l - rSTCY 


AP-5061-R 
Appendix  C 

38 


Kote  that  to  make  the  terms  a  function  of  forward  valoeltir 
(to  die  out  at  high  speed)  the  terni  (1  +  e  -  )  li 

Qr 

added  to  the  above  terma*. 

This  takes  Into  account  the  time  response  of  the  water 
surface.  At  high  speed  the  GEM  base  Is  above  a  given  point  on 
the  water  for  such  a  short  time  that  the  surface  may  deflect  only 
a  small  amount  (compared  with  Its  static  deflectidh)  before  the  GEM 
has  passed  by. 
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ANALOG  COMPUTER  DIAGRAM 
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